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Abstract 

The coupling of the metric to an incoherent dust introduces into spacetime 
a privileged dynamical reference frame and time foliation. The comoving 



o . 

coordinates of the dust particles and the proper time along the dust world- 

G\ 
O 

t^J- \ Hamiltonian constraint can be resolved with respect to the momentum that 

is canonically conjugate to the dust time. Imposition of the resolved con- 
straint as an operator restriction on the quantum states yields a functional 
i , Schrodinger equation. The ensuing Hamiltonian density has an extraordi- 

■ 

(3JT), nary feature: it depends only on the geometric variables, not on the dust 

coordinates or time. This has three important consequences. First, the func- 
tional Schrodinger equation can be solved by separating the dust time from 
j_j ■ the geometric variables. Second, the Hamiltonian densities strongly commute 

and therefore can be simultaneously defined by spectral analysis. Third, the 
standard constraint system of vacuum gravity is cast into a form in which 
it generates a true Lie algebra. The particles of dust introduce into space a 
privileged system of coordinates that allows the supermomentum constraint 
to be solved explicitly. The Schrodinger equation yields a conserved inner 
product that can be written in terms of either the instantaneous state func- 
tional or the solutions of constraints. Gravitational observables admit a sim- 
ilar dual representation. Examples of observables are given, though neither 
the intrinsic metric nor the extrinsic curvature are observables. Disregarding 
the standard factor-ordering difficulties, the introduction of dust provides a 
satisfactory phenomenological approach to the problem of time in canonical 
quantum gravity. 
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I. INTRODUCTION 



The Dirac constraint quantization of vacuum Einstein gravity yields the Wheeler-DeWitt 
equation for the quantum state of the intrinsic three geometry of space . One can view 
this equation as a statement that only two out of three independent components of the 
intrinsic geometry are dynamical. The third component is an intrinsic time that specifies 
the location of space as a hypersurface in spacetime. The Wheeler-DeWitt equation is then 
interpreted as an evolution equation for the state in the intrinsic time. 

The Wheeler-DeWitt equation is a second-order variational differential equation. The 
space of its solutions carries no obvious Hilbert space structure [@,|J. This has prompted 
numerous attempts aimed at replacing the Wheeler-DeWitt equation by a first-order 
Schrodinger equation. In order to do that, one should identify the intrinsic time at the 
classical level, solve the Hamiltonian constraint for the momentum conjugate to time, and 
impose the resolved constraint as an operator restriction on the quantum states. Unfortu- 
nately, there is no natural candidate for the intrinsic time, and the procedure is beset by a 
number of conceptual and technical difficulties [f| . 

Intrinsic clocks are strange contraptions. From the early days of general relativity, most 
researchers felt that spacetime intervals are to be measured not by internal clocks but by 
material devices, analogous to the measurement of electromagnetic fields by test charges. 
At the most basic level, matter splits the spacetime manifold into space and time. The 
phenomenological medium introduced for this purpose is the reference fluid The par- 

ticles of the reference fluid identify the points of space, and clocks carried by these particles 
identify the instants of time. In this way, the fluid fixes the reference frame (the space) and 
the time foliation (the time). In that frame and on that foliation, the entire intrinsic metric, 
not just two selected components of the intrinsic geometry, becomes dynamical. 

The reference fluid is traditionally considered to be a tenuous material medium whose 
back reaction on the geometry can be neglected. There are just enough fluid particles to 
discern the points of space from one another, but not enough to disturb the geometry. 
Instead of deriving the motion of the fluid from its action, one encodes it in coordinate 
conditions. Unfortunately, such a standpoint makes it difficult to view the reference fluid as 
physical matter. 

There are two alternative routes through which the reference fluid can be turned into 
a physical system. The first is to impose the coordinate conditions before variation by 



adjoining them to the Hilbert action with Lagrange multipliers [ 1QU11 ■ Equivalently, a 



Lagrangian might be devised to enforce these conditions [|12|,[nj. The additional terms 
in the action can be interpreted as a matter source coupled to the gravitational field. A 
disadvantage of this approach is that simple coordinate conditions often lead to contrived 
matter sources. The second route is to select a realistic material medium for the fluid and 
to describe its properties by a physical Lagrangian P^,|2| . Again, by adding this Lagrangian 
to the Hilbert Lagrangian, the fluid becomes coupled to gravity. The challenge is to find 
a medium whose space and time variables can be turned into canonical coordinates that 
lead, upon the Dirac constraint quantization, to a simple Schrodinger equation for canonical 
quantum gravity. 

Our goal is to show that the most simple type of matter conceivable, incoherent dust, 
is just what one needs to obtain a satisfactory framework for interpreting quantum gravity. 
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The particles of dust follow timelike geodesies, and in this sense dust is quite special: Dust 
is the 'purest' form of matter, matter that interacts with itself only via the gravitational 
field. 

The canonical scheme offered by incoherent dust has an extraordinary feature: The 
many-fingered Hamiltonian of the ensuing functional Schrodinger equation does not depend 
on the dust variables. This fact has a number of important consequences. First, the Hamil- 
tonian densities strongly commute (modulo factor-ordering difficulties), and thus can be 
simultaneously defined by spectral analysis. Second, the functional Schrodinger equation 
can be solved by separating the dust time from the gravitational variables. Third, the fact 
that the Hamiltonian depends only on the gravitational variables leads to the discovery of 
a form for the constraint system of vacuum gravity in which the constraints generate a true 
Lie algebra. 

The dust particles also endow space with a privileged, physical system of coordinates. 
Canonical variables such as the spatial metric can be specified as tensors in these dust 
coordinates, and in this way the combined system of gravity and dust can be expressed in a 
form that does not rely on arbitrary systems of spatial coordinates |T5| . As a consequence, 
the supermomentum constraint, which is the canonical generator of spatial diffeomorphisms, 
becomes trivial and is exactly solvable in the quantum theory. 

All of the above features combine into an intriguing scheme for interpreting canonical 
quantum gravity. Although we focus on the case of gravity coupled only to dust, the scheme 
can be applied in the presence of other matter couplings as well — as long as dust is included 
in the list of matter fields, all of our basic conclusions remain intact. Here is how the scheme 
works: 

We design a spacetime action whose variation yields the Euler equations of motion for 
an incoherent dust (Section 2a). All of the state variables have a direct physical significance. 
In particular, four of these variables are the proper time and the dust frame (Section 2b); 
these are the desired tools for interpreting quantum gravity. The action is unchanged by 
a subgroup of diffeomorphisms in the dust spacetime manifold (Appendix C). These global 
symmetries yield an infinite set of Noether charges (Section 2e) that are conserved by virtue 
of continuity equations (Section 2c). By varying the dust action with respect to the spacetime 
metric, the standard energy-momentum tensor is obtained (Section 2d). The dust action 
can be cast into several equivalent forms (Appendix B). 

The ADM decomposition of the dust action, followed by a Legendre dual transformation, 
yields the canonical description of dust (Section 3a). The canonical coordinates coincide with 
the dust time and frame variables, while the conjugate momenta are simply related to the 
mass density and velocity of the dust. The canonical form of the dust action can be derived 
from the action for a continuously infinite number of relativistic point particles (Appendix 
D), and is a special case of the canonical action for arbitrary perfect fluids jnj. The Noether 
charges are represented by linear functionals of the dust momenta (Section 3b). 

Dust is coupled to gravity by adding its action to the Dirac-ADM action. The dynamics 
is then generated by the standard super-Hamiltonian and supermomentum constraints on 
the phase space of the coupled system (Section 4a). The super-Hamiltonian and super- 
momentum constraints can be brought to a form in which they are resolved in the dust 
momenta. These new constraints H^(x) and H^k{x) generate an Abelian algebra (Section 
4b) that allows us to represent diffeomorphisms in the dust spacetime by canonical trans- 
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formations (Section 4c). 

The new super-Hamiltonian constraint H^(x) generates dynamics along the flow lines 
of the dust, while the new supermomentum constraint H^ k {x) generates motion along the 
hypersurfaces of constant dust proper time. The true Hamiltonian density in the super- 
Hamiltonian constraint is the square root of a quadratic combination G(x) of the grav- 
itational super-Hamiltonian and supermomentum. As such, it does not depend on the 
dust variables. The expressions G(x) have vanishing Poisson brackets among themselves 
(Appendix E). When complemented by the gravitational supermomenta, they provide an 
alternative set of constraints for vacuum gravity which generate a true Lie algebra (Section 
4d). 

The dust frame variables are fields Z k (x), k = 1,2, 3, whose values at a given space point 
x constitute a set of Lagrangian coordinate labels for the dust flow line that passes through 
x. Z can be viewed as a mapping from the space manifold to the 'dust space' manifold whose 
points z are the individual flow lines (Section 2b). Correspondingly, Z induces mappings 
of tensor fields. In particular, the spatial metric, the dust proper time variable, and their 
canonical conjugates can be mapped from space to the dust space. This mapping of tensors 
is completed into a canonical transformation by replacing the momentum conjugate to the 
dust frame Z k [x) with the new supermomentum generator H^(x) (Section 5a). Then the 
constraint H^ k (x) = is simply the statement that the momentum conjugate to Z k [x) 
should vanish. Because the new super-Hamiltonian constraint (x) is independent of the 
dust frame Z and its conjugate, it can be expressed in terms of the new canonical variables by 
mapping the space scalar density H^(x) to the dust space. This turns the super-Hamiltonian 
constraint into a dust space scalar density H^(z). The true Hamiltonian density in this 
constraint is the square root of a dust space density G(z) that is constructed solely from 
the gravitational variables. The Noether charges also can be expressed as functionals of the 
new canonical variables (Section 5b). 

We quantize gravity coupled to dust by imposing the supermomentum constraint 
Hffc(x) = and the new super-Hamiltonian constraint H^(z) = as operator restrictions 
on the states (Section 6). The supermomentum operator constraint implies that the state 
functional is independent of Z, while the new Hamiltonian constraint yields a functional 
Schrodinger equation in the dust time (Section 6a). Since the Hamiltonian density in this 
equation does not depend on the dust variables, the dust time can be separated from the 
gravitational variables. This leads to a time-independent functional Schrodinger equation 
(Section 6b). If the operators G(z) can be factor ordered and regularized so that they com- 
mute, then they possess joint eigenfunctionals (Section 6c). We show that the spectrum of 
G(z) is continuous and degenerate (Appendix F). The square-root Hamiltonian density can 
be defined by spectral analysis on the Hilbert space 7i + spanned by the eigenfunctionals la- 
beled by positive eigenvalues G(z) > 0. The Noether charges appear in the quantum theory 
as generators for unitary operators that transform the quantum states by the corresponding 
symmetries (Section 6d). The quantum theory of gravity and dust can be constructed using 
either the original set or the new set of canonical variables (Section 6g). 

The Schrodinger equation yields a conserved inner product that can be written either 
in terms of the instantaneous state functionals, or as a product in the space of solutions 
(Section 6e). Gravitational observables admit a similar dual representation (Section 6f). 
Observables must keep the state in the Hilbert space H + . Examples of such observables are 
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given; unfortunately, neither the intrinsic metric nor the extrinsic curvature are gravitational 
observables. A similar difficulty occurs already in a simple model system, namely, the 
relativistic particle in the proper-time formalism (Appendix G). 

The classical Noether charges generate canonical transformations that transform the 
system from one dust time to another (Section 3b). They carry a solution of the constraints 
and of the Hamilton equations of motion again into a solution. If the dust does not rotate, 
the solution can be transformed to the time foliation that is orthogonal to the dust frame 
(Section 7a). The restriction to non-rotating dust in the frame-orthogonal time foliation can 
be enforced by an additional first-class constraint. This constraint can be used to eliminate 
the inconvenient square root in the Hamiltonian density H^(x) (Section 7b). The resulting 
Hamiltonian constraint H^o(x) underlies our previous study of Gaussian reference fluids 
||10|| . Due to factor-ordering ambiguities, the quantum theories constructed from H^(x) and 
Hi (x) do not necessarily coincide (Section 7c). 

Our notation is summarized in Appendix A. 

The introduction of dust alleviates several notorious problems of time in vacuum gravity: 

• When passing to the Schrodinger equation, one can rely on the privileged time and space 
brought into spacetime by the state of a material system. 

• The dust time is a spacetime scalar. Its value does not depend on the hypersurface passing 
through a given event. No intrinsic time in vacuum gravity has this property. 

• The Hamiltonian densities have vanishing Poisson brackets. If — and this is an overwhelm- 
ing if — a factor ordering and regularization can be found for which the corresponding op- 
erators G{z) commute, then the Hamiltonian density operators can be defined by spectral 
analysis and the evolution of states will be independent of the foliation connecting the final 
and initial embeddings. One important problem remains to be solved: What operators, if 
any, correspond to the measurements of the intrinsic metric and extrinsic curvature on an 
embedding specified with respect to the dust time and the dust frame? 

To the extent limited by our cautionary remarks, the introduction of incoherent dust 
provides a satisfactory phenomenological approach to the problems of time in canonical 
quantum gravity. The corresponding problems at the level of vacuum gravity remain, of 
course, unresolved. 

Preliminary accounts of our work have been given in Refs. [|l7jj . 

II. INCOHERENT DUST 

A. Spacetime Action and the Euler Equations 

Dust is described by eight spacetime scalars, T, Z k , M, and Wk- (Indices k, £ from the 
middle of the Latin alphabet take the values 1, 2, 3.) The physical interpretation of these 
state variables follows from the analysis of the action and the ensuing equations of motion. 
In order to describe the structure of the action in physical terms, we shall anticipate the 
results of this analysis. The values of the variables Z k are the comoving coordinates of the 
dust particles, and the value of the variable T is the proper time measured along the particle 
flow lines. The scalars Z K = (T, Z k ) are assumed to be four independent functions of the 
spacetime coordinates y a , a = 0, 1, 2, 3: 
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det (Z K :(X ) ^ . (2.1) 

The four-velocity U a of the dust is defined by its decomposition in the cobasis Z K a as 

U a = -T a + W k Z\ a . (2.2) 

The three spatial components W k of the four-velocity in the dust frame Z k are also state 
variables. Equation (2.2) expresses the four-velocity one-form U = U a dy a as a Pfaff form 



U = -dT + W k dZ k (2.3) 

of the seven scalar fields T,Z k and W k . The eighth state variable, M, is the rest mass 
density of the dust. 

The dust action is a functional of the dust variables and of the spacetime metric j a /3, 

S D [T,Z k ;M,W k ;^]= [ d 4 yL D (y), (2.4) 

J M 

with the Lagrangian density 

L D = ~\j\ 1/2 M ^XJJJp + l) . (2.5) 

It can be cast into many alternative forms (Appendix B). The matter equations of motion 
are obtained by varying the dust action with respect to the state variables M, W k , T, and 
Z k : 







5S 



D 
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\\l\ 1/2 {l aP U a U, + l) , (2.6) 



= — - = -\ 1 \ 1/2 MZ k , a U a , (2.7) 
dW k 

° = -§r = -W A,0 °)*' (2 ' 8) 

£ CD 

= j^ = (\l\ 1/2 MW k U^ a . (2.9) 
These equations corroborate our interpretation of the dust variables: 



B. Dust As a Standard of Space and Time 

Equation (2.6) ensures that the four-velocity U a is a unit timelike vector field. According 
to Eq. (2.7), the three scalar fields Z k are constant along the flow lines U a and therefore 
their values z k can be interpreted as Lagrangian coordinates for the dust. That is, the values 
z k = Z k (y) of the mappings Z k : Ai — > M serve as labels specifying which flow line passes 
through a given spacetime event y. Such a set of Lagrangian coordinates can be generated 
by choosing an arbitrary spacelike hypersurface £ and specifying a coordinate system z k on 
that hypersurface. Then each flow line is labeled by the coordinates of the point where it 
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intersects the hypersurface. Inversely, knowing the congruence S of flow lines, we can define 
the mappings Z k of Ai into M by assigning to each event y in Ai the coordinates z k of the 
point in £ where the flow line that passes through y intersects E. 

The congruence S of dust flow lines introduces a privileged reference frame into the 
spacetime manifold Ai. This congruence can be viewed as an abstract three-dimensional 
space, the 'dust space', whose points are the individual flow lines. As a manifold, the dust 
space S is isomorphic to any spacelike hypersurface. (We assume that the spacetime (Ai, 
7) is globally hyperbolic with Ai = M x £.) A system of coordinates on that hypersurface 
induces a system of coordinates z k on S. It may be impossible to cover the dust space with a 
single coordinate chart, so the Lagrangian coordinates z k are generally defined only in open 
subsets of the dust space. Correspondingly, the scalar fields Z k are generally defined only 
in open tubes of the space-time manifold. Our interpretation of the fields Z k presupposes 
that their values constitute a good set of Lagrangian coordinates, i.e., that within the 
appropriate open tubes each flow line carries a unique set of labels z k . As a consequence, 
the three gradients Z k ^ a must be linearly independent covectors. This is ensured by the 



assumption ( |2.1| ) 



Our discussion reveals that the three scalar fields Z k (y) are not necessarily globally 
defined on the whole spacetime manifold Ai. However, these fields appear in the action (2.4) 
only in the combination WkZ k ^ a , so this defect can be overcome by a slight reinterpretation 
of the formalism. Let Z : Ai — > S denote a globally defined mapping of the spacetime 
Ai into the dust space S, and view Z k as the composition of the mapping Z with a local 
coordinate chart mapping on S. With the indices k and a interpreted as abstract tensor 
indices, Z >a is then globally defined as the mapping Z k )CC : At — > TS ®T*Ai. Likewise, Wk 
is defined as the mapping Wk '■ Ai — > T*S where, again, k is viewed as an abstract index. 
In this way, WkZ k tOC is understood to be a globally defined element of T*Ai. 

The dust enables us to introduce not only the privileged matter space, but also a matter 
time. By multiplying Eq. (2.2) by U a and using the field equations (2.6)-(2.7), we learn 
that 

£uT = T a U a = 1 , (2.10) 

i.e., that T is the proper time between a fiducial hypersurface T = and an arbitrary 
hypersurface T = const along the flow lines of the dust particles. The fiducial hypersurface 
can be chosen to be spacelike. The foliation of the spacetime manifold Ai by hypersurfaces 
of constant T can be viewed as an abstract one-dimensional space, the 'dust time' T , 
whose points (instants) are the individual hypersurfaces. As a manifold, the dust time T is 
isomorphic to M. The value r = T(y) of the time mapping T : Ai — > T can serve as the 
time coordinate of an instant of time that passes through the event y. 

In our analysis, we must assume that the three covectors Z >a are linearly independent. 
That the fourth covector, T Q , is linearly independent of Z k >a follows from the equations 
of motion (2.7) and (2.10). As a result, (T a ,Z k >a ) form a basis in T*M, Eq. Q). The 
four-velocity U a of the dust is a privileged timelike future-pointing vector. Equation (2.10) 
tells us that the time function T(y) grows from the past to the future. The covectors Z k >a 
are orthogonal to the four-velocity U a of the dust and hence spacelike. These relations are 
sketched in Fig. 1. 
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Because (T a , Z k ^ a ) are linearly independent, the values r and z k of the four mappings T 
and Z k can be used as coordinates on Ai. By inverting the mappings 

r = T(y), z = Z(y) (2.11) 

we obtain the mapping 

T-.TxS^M by (r,z) ^y = T(t,z) . (2.12) 

The four vectors 

U a := T Q r and T a fe (2.13) 

form a basis in TM. dual to the cobasis 

Z K ta = (T a ,Z k ia ) (2.14) 

in T*M.. The basis (2.13) and the cobasis (2.14) satisfy the standard orthonormality and 
completeness relations. As a consequence of the equations of motion (2.7) and (2.10), these 
basis and cobasis vectors are conserved along the flow lines of U a \ that is, £\jT a = = 
£\jZ k jQ and £uU a = = £\jY a )k . The basis vector U a := T° T coincides with the four- 
velocity of the dust and the remaining three vectors, T"^, are tangent to the T = const 
hypersurfaces. These relations are again sketched in Fig. 1. 

To summarize, any solution of the dust equations fl2.6|) and (2.7) which satisfies the 
condition fl2.1p describes a motion of the dust that allows the spacetime manifold A4 to be 
split into the space and time manifolds, S and T. An instant of T and a flow line of S 
intersect at a unique event y G M. The spacetime M. is thus a Cartesian product T x S of 
the space and time manifolds. The mapping T(y) x Z(y) takes us from M. to T x S, the 
inverse mapping T(r, z) takes us from T x S to M.. In this sense, dust serves as a standard 
of space and time. 



C. Conservation Laws 

Let us now turn to the second set of Euler equations. Equation (2.8) is a continuity 
equation for the rest mass current 

J Q := MU Q . (2.15) 

The total mass is thus conserved from one spacelike hypersurface to another. Similarly, the 
Euler equation (2.9) is a continuity equation for 

J fc Q := W k J a , (2.16) 

which we refer to as the momentum current. These two continuity equations imply 

£uW k =W Ka U a = , (2.17) 

so that Wk stays the same along each flow line. As a result, the velocity W k can depend on 
y only through the Lagrangian coordinates z = Z(y): 

W k (y) = W k (Z(y)) . (2.18) 

This shows that W k can be viewed as a covector, or one-form, in the dust space S. 
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D. The Energy Momentum Tensor 



Dust is coupled to the metric field by adding its action S D to the Hilbert action S , 
which is given by 

S G h a p]= I d 4 y L G , L G = \ 1 \ 1 l 2 R{y ]laP ] . (2.19) 

J M 

The variation of the total action S with respect to 7 a/ g yields the Einstein law of gravitation 

R a(3 - -i? 7 Q/3 = -T af3 (2.20) 
with the energy-momentum tensor 

T *t3 := 2 | 7 |i/2 = MU a UP + - 1 aP M( 1 ^UJJ v + 1) 

Ola/3 2 

« MU a U p . (2.21) 

Modulo the normalization equation (|2.6| ), this tensor has the form required for the description 
of incoherent dust with rest mass density M and four- velocity U a . 

By the well-known argument, the invariance of the action (2.4)-(2.5) under spacetime 
diffeomorphisms (p, namely, 

S D [tp*T, <p*Z k ; ip*M, ip*W k - V * laf5 } = S D \T,Z k -M,W k - la!i \ , (2.22) 

leads to the identity 

M"*V^ , S -^T a + ^ + + . (2.23, 

It follows that if the dust satisfies the field equations (2.6)-(2.9), the energy-momentum 
tensor is covariantly conserved: 

V^Tf = . (2.24) 

The conservation law (2.24) also follows independently of the Euler equations (2.6)-(2.9) 
from the Einstein law of gravitation (2.20) through the contracted Bianchi identities. The 
identity (2.23) shows that if we postulate the conservation law (2.24), the Euler equations 

(2.6) -(2.9) are not independent. It suffices to impose the first set of equations, (2.6) and 

(2.7) , obtained by varying the multipliers M and W k - The second set of Euler equations, 

(2.8) and (2.9), obtained by varying the dynamical variables T and Z fc , then follows because 
(T a , Z k >a ) is a regular cobasis. 

The first and second sets of Euler equations have a complementary character. The first 
set, obtained by varying the multipliers M and W k , tells us that the dynamical variables 
T and Z k have the properties (2.10) and (2.7) of the proper time and the Lagrangian 
coordinates of the dust. The second set, obtained by varying the dynamical variables T 
and tells us that the multipliers M and W k satisfy conservation laws (2.8) and (2.9) (or 
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(2.17)). In the canonical formalism, the state variables T and Z k are turned into canonical 
coordinates, while the Lagrange multipliers M and are combined to form the conjugate 
momenta. 

From the energy-momentum conservation (2.24) and the mass conservation (2.8) it fol- 
lows that the flow lines of the dust particles are geodesies, 

U p V p U a = . (2.25) 

This demonstrates that the dust particles do not act on one another by mechanical forces. 

E. Symmetries and Noether Charges 

The action (2.4)-(2.5) depends on the variables T, Z k , and Wk only through the combi- 
nation that occurs in the Pfaff form (2.3). It follows that invertible transformations of these 
variables that leave the Pfaff form unchanged constitute symmetries of the dust action. In 
Appendix C we identify two particular types of symmetry transformations that mix the dust 
spacetime variables among themselves. The first is expressed as 

Z k = s *( Z m') ) Wk , = 3^(Z m ') W e . (2.26) 

These transformations are diffeomorphisms of the dust space S, where the velocity Wk 
transforms as a covector. Under the transformations (2.26), the dust worldline that passes 
through the spacetime point y has its Lagrangian coordinate labels changed from z k = Z k (y) 
to z k ' = Z k \y). 

The second type of symmetry transformation is expressed as 

T' = T + &(Z e ) , W v = W k + &, k (Z e ) . (2.27) 

The dust time T is not entirely determined by the dust flow lines since the initial hypersurface 
can be chosen arbitrarily. Once an initial hypersurface is chosen, all other hyper surf aces of 
the time foliation are determined, according to Eq. (2.10), by Lie propagation of the initial 
hypersurface along the flow lines of U a . If the initial hypersurface is changed, the whole 
foliation T is changed. Equation (2.27) connects the time function T(y) of the old foliation 
T(y) = t = const with the time function T'(y) of the new foliation T'(y) = r' = const. It 
also determines how the projection Wk of U a onto the leaves of the time foliation is affected 
by a different tilt of the leaves with respect to the flow lines. 

The transformations (2.26)-(2.27) are global symmetries of the action, because they do 
not depend explicitly on y and therefore cannot be made to vanish in a compact region of 
M.. (Gauge symmetries, on the other hand, vary with y.) These symmetry transformations 
carry one solution of the field equations into another solution. Indeed, it is easy to check 
that if the field equations ( |2.7| ) and ( |2.9| ) — the only ones that contain the relevant state 
variables — are satisfied for the original variables, they are also satisfied for the transformed 
variables. 

According to Noether's theorem, to each one-parameter group of symmetry transforma- 
tions connected to the identity there corresponds a conserved charge. Let ip r (y) denote a 
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collection of field variables and ip' r = \l/ r (V ; , cr) denote a one-parameter group of symmetries 
with the identity at a = 0. The corresponding conserved charge Q is given by the integral 



of the Noether current 



J a : = 



dL dV r (il>,a) 



da 



(2.28) 



(2.29) 



a=0 



The integral in Eq. (2.28) is over an arbitrary spacelike hypersurface Y(E) specified by an 
embedding 



Y : £ ^ M by x eZ^y = Y(x) 
and carrying the oriented surface element 

rfS a := 5 a p lS Y^ h Y\ c Y\ d dx b A dx c A dx d 



(2.30) 



(2.31) 



(Note that J a is a spacetime vector density of weight 1, and dE a is a covector density of 
the opposite weight —1 due to the presence of the alternating symbol S^ a /3yS-) 

The symmetry transformations (2.26)-(2.27) contain arbitrary functions S k (Z e ') and 
@(Z e ). The Lie algebra of their one-parameter subgroups is thus labeled by the arbitrary 
functions 



dS k (Z e ',a) 



da 



&(Z e ) :-- 



de(Z e ,a) 



CT=0 



da 



(2.32) 



(7=0 



For the symmetry transformations (2.26), the corresponding Noether currents (2.29) are 

J a [£] = \l\ 1/2 J^\Z) , (2.33) 

where J% is the momentum current (2.16). For the symmetry transformations (2.27), the 
corresponding Noether currents are 



J«[0] = | 7 |V2j^(Z) , 



(2.34) 



where J a is the mass current (2.15). By virtue of the field equations (2.7)-(2.9), the Noether 
currents satisfy the continuity equations 



J%[£] = o = r,M ■ 

The Noether charges (2.28) can be written in the form 



Q[t] 



d*x \g\ 1 / 2 n a JZ£ k (Z(x)) , 



Q[#] = - [ d 3 x \g\^ 2 n a J a ^(Z(x)) , 



(2.35) 

(2.36) 
(2.37) 
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where \g\ l l 2 d?x is the proper volume element of Y(H) and n a is the future-pointing unit 
normal. 

In order to interpret the charge (2.37), choose "d to be unity for all Z k in some ball B 
in the dust space S, and zero outside B. Then Q\&\ becomes a proper volume integral over 
the subspace of the hypersurface Y(E) that contains flow lines in £>, with the integrand 
M(—n a U a ). The factor —n a U a is just the relativistic 'gamma factor' characterizing a boost 
from the Lagrangian observers whose four-velocities are U a to the Eulerian observers on 
Y{T) whose four- velocities are n a . Thus, Q{&\ is the rest mass in B as seen by the Eulerian 
observers. Conservation of the charge (2.37) expresses the conservation of relativistic mass 
within a flow tube defined by the bundle of flow lines crossing the ball B. For B — E, 
conservation of the total relativistic mass is obtained. 

Similarly, choose $, k = 5 k in the ball B, and zero outside B. The charge (2.36) has the 
meaning of minus the component of dust momentum along the coordinate line z l in the 
T = const surface within B. 



III. CANONICAL DESCRIPTION OF DUST 

A. The Canonical Action 

The ADM decomposition of the dust action follows the familiar pattern. Foliate space- 
time into spacelike hypersurfaces by 

Y-.lRxZ^M by (t, x) i-> y = Y(t, x) , (3.1) 

a one-parameter family of embeddings of space £ into spacetime M.. In local coordinates 
x a on E and y a on M, the foliation (3.1) is represented by 

{t,x a ) ^y a = Y a {t,x a ) . (3.2) 

A transition from one leaf of the foliation to another is described by the deformation vector 
N a := Y a . Its decomposition 

N a := Y a = N L n a + N a Y a A (3.3) 

into the unit normal n a and tangential Y a >a directions yields the lapse function N 1 - and the 
shift vector iV a . The spacetime metric 7 a/ 3 is determined by 

rf& = - n * n P + g ab Y a , a Y^ b , M 1 / 2 = N^ 1 ' 2 , (3.4) 

where g ab is the inverse of the spatial metric g a b induced on E. 

Scalar fields on A4, such as Z k , can be pulled back to M x E by the mapping (3.1). 
Then, according to Eq. (3.3), 

Z\ a n a = {N L )-\Z k - Z k A N a ) . (3.5) 

Similar equations hold for the scalar fields T, M, and Wk- From a global point of view, 
the mapping (3.1) is used to pull back the mapping Z : M. — > S to 1R x E, which yields a 
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one-parameter family of mappings Z t : S — > 5 from space E into the dust space 5. (We will 
often suppress the t-dependence and simply write Z : E — > S.) With the indices and a 
viewed as abstract indices, Z k denotes a one-parameter family of mappings from E to TS. 
Likewise, Z k ^ a denotes a one-parameter family of mappings from E to TS <g> T*E and W k 
denotes a one-parameter family of mappings from E to T*S. 

The results above enable the dust action (2.4)-(2.5) to be written as an integral over 
R x E, 

S D [T,Z k ; M,W k ; g^N^^N*] = f dt f d 3 x L D , (3.6) 



where 



° l \g\ l ' 2 M Qzi-T + W k Z k -U a N a ) 2 - N\g ab U a U b + l)) , (3.7) 

is the dust Lagrangian density and U a := U a Y a A = —T a + W k Z k A . 
Let us introduce the momentum P conjugate to T by 

dL D 

P := — ^ = |0| 1 / 2 M(JV- L )- 1 (T - W k Z k + U a N a ) . (3.8) 
dT 

The physical meaning of P is revealed when we identify the relativistic 'gamma factor' 
—U a n a characterizing a boost from the Lagrangian observers whose four-velocities are U a , 
to the Eulerian observers on Y(E) whose four-velocities are n a . From Eq. (2.2) for U a and 
Eqs. (3.5) for Z k ^ a n a and T Q n a we obtain 

- U a n a = (N^- 1 (f - W k Z k + U a N a ) . (3.9) 

This shows that the canonical momentum P is equal to the projection 

P=\g\ 1/2 (-J a n a ) (3.10) 

of the rest mass four-current (2.15) into the four-velocity of the Eulerian observers. There- 
fore P is the density of rest mass on Y(E), per unit coordinate cell d 3 x (rather than per unit 
proper volume dV = \g\ 1 ^ 2 d 3 x), as measured by the Eulerian observers. Equation (3.10) also 
shows that the momentum P is necessarily positive, 

P>0, (3.11) 

since J a and n a are both future-pointing timelike vectors. 

Perform now the Legendre dual transformation (T, T) — > (T, P) while leaving the vari- 
ables Z k , M, and W k untouched. At the same time, introduce the abbreviation 

P k := -PW k (3.12) 

for the Lagrange multiplier W k scaled by —P. The resulting Hamiltonian PT — L D is 

pf - L D = -P k Z k + N L H I l + N a H° , (3.13) 
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with 



H° := -PU a = PT a + P k Z\ a (3.14) 

and 



Equation (3.13) reveals that the Legendre transformed action 

S D [T,P; Z k ,P k ; M; ^, JV" 1 , JV a ] 

= [ dt [ d 3 x(PT + P k Z k - N L Hl - N a H°) (3.16) 

J IFt J S 

is in canonical form not only in the variables T and P, but also in the variables Z k and P k . 

The variable (3.12) is thus the momentum canonically conjugate to Z k . As a pendant to 
Eq. (3.10), the momentum P k of Eq. (3.12) can be expressed as 

P k = -\g\ 1/2 {-J%n a ) , (3.17) 

and interpreted as (minus) the momentum density along T a >k per unit coordinate cell d 3 x, 
measured by the Eulerian observers on Y(E). 

By varying the action (3.16) with respect to M we get 

^7 = = M = , * , = (3.18) 
5M \9\ 1/2 JP*+g«>>H?H b D 



The positive square root is dictated by the positivity of M. Substituting the expression 
(3.18) back into iff gives 



and eliminates the mass multiplier M from the action. Equations (3.12) and (3.18) connect 
the multiplier variables M, W k of the spacetime action with the canonical momenta P, P k - 
The spatial densities and are, respectively, minus the momentum density and 
energy density (in the coordinate basis dx a ) of dust as seen by the Eulerian observers n a on 
Y(E). The canonical form (3.14) of is dictated by its role as the canonical generator 
of spatial diffeomorphisms for the spatial scalars T, Z k and their conjugates P, P k . Very 
few field theoretical systems have a square-root energy density like iff. (Two such systems 
known to us are the Born-Infeld electrodynamics |19[ and the extrinsic time reference fluid 



The square-root structure of the dust Hamiltonian (3.19) arises from the square-root 
structure of the Hamiltonian of a single relativistic particle. This connection is demonstrated 
in Appendix D, which makes use of results from Sec. 5. 

Because the dust action contains no derivatives of the metric, the ADM decomposition 
of the total action leads to the super-Hamiltonian and supermomentum constraints 
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H± := Hf + Hi = , (3.20) 

H a :=HC + H° = 0, (3.21) 

whose gravitational contributions Hf and H£ have the same form as in vacuum general 
relativity, namely, 

Hl(x; g ab ,p ab ] = G abcd (x; 9) p a \x)p cd (x) - \g\ l ' 2 R{x- g] , (3.22a) 

Gabcd = ^\g\~ ll2 {gac9bd + 9ad9bc ~ 9ab9cd) , (3.22b) 

H°(x;g ab ,p ab ] = -2D b p b a (x) . (3.23) 



The total constraints (3.20) and (3.21) satisfy the usual Dirac 'algebra' (Sec. 4a). 

B. Symmetries and Charges 

By using Eqs. (3.10) and (3.17), the Noether charges (2.36) and (2.37) can be expressed 
as functionals of the canonical variables, 

Q[£\ = ~ I d s x e(Z(x))P k (x) , (3.24) 
QVd] = [ d 3 x &(Z(x))P(x) . (3.25) 

The global symmetries of Section 2e appear in the canonical formalism as canonical transfor- 
mations of the dynamical variables generated through the Poisson brackets by the Noether 
charges (3.24) and (3.25). For example, by applying the canonical transformation generated 
by the charge (3.25) to the dynamical variables T(x) and 

W k (x) := -P k (x)/P(x) , (3.26) 

the infinitesimal version of Eq. (2.27) is recovered: 

{T(x),Q[ti]} = ti(Z(x)) , {W k (x),Q[0]} = #k( z (x)) ■ ( 3 - 27 ) 

Likewise, the canonical transformation generated by the charge (3.24) yields the infinitesimal 
version of Eq. (2.26). 

The phase space functionals (3.24), (3.25) have vanishing Poisson brackets with the 
constraints (3.20) and (3.21), so they indeed generate symmetries of the theory. By the 
same token, they are constants of motion. The symmetry generators (3.24) and (3.25) close 
under the Poisson brackets according to the equations 

{Q[£iW[£ 2 ]} = a^] , (3.28a) 

{Q[&Qm = Q[£i*] , (3.28b) 

{Q[0i],Q[0 2 ]} = O. (3.28c) 

Thus, the Q's form a realization of the Lie algebra of the subgroup (2.26)-(2.27) of diffeo- 
morphisms of the dust spacetime T x S. 
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IV. THE CONSTRAINT ALGEBRA 



A. The Dirac Algebra of Old Constraints 

The constraints (3.20)-(3.21) of dust coupled to gravity satisfy the usual Dirac 'algebra' 

{H ± (x),H ± (x')} = g ab (x)H b (x)5 >a (x,x') - (x <-> x') , (4.1) 
{H ± (x),H a (x')} = -H ± (x')5 >a ,(x,x') , (4.2) 
{H a (x), Hf,(x')} = Hb(x)5 ja (x,x') — (ax «-> bx') . (4.3) 

Equation (4.3) tells us that H a (x) represent the generators of the spatial diffeomorphism 
group DiffS. Equation (4.2) expresses the fact that the super-Hamiltonian H±(x) is a scalar 
density of weight 1 under DiffS. The right-hand-side of Eq. (4.1) contains the dynamical 



variable g a b(x). This implies that the total system (4.1)-(4.3) is not a true Lie algebra [20 

The gravitational super-Hamiltonian H^(x) and supermomentum H^(x) close according 
to the same relations (4.1)-(4.3) as do the complete constraints H±(x) and H a (x). From 
here, one can deduce the Poisson brackets among the matter parts H±(x) and H®(x) of 
the constraints [21|. The momentum densities H®(x) still represent DiffS, Eq. (4.3). For 



systems with non-derivative gravitational coupling (which is the case for dust), the energy 
density H®(x) depends only on the metric g a b(x), but not on its derivatives. It follows 
that H±(x) and H^(x') close in the same manner, Eq. (4.1), as the total constraints H±(x) 
and H±(x'), or as the gravitational expressions H^(x) and H^_(x'). The Poisson bracket 
{iff (x), (%')} does not follow the pattern (4.2) because iff (a/) does not act on the 
metric in iff (x). The concrete form of this bracket will not concern us here. 
The energy density iff (a?) of the dust is the square root of the expression 

D(x) := P\x) + g ab (x)H°(x)H b D (x) . (4.4) 

From the closing relations (4.1) for H D (x) it follows that 

{D(x), D(x')} = AD(x)H Da (x)5 ta (x, x') - (x «-> x') . (4.5) 



B. The Abelian Algebra of New Constraints 

The fact that the Dirac 'algebra' is not a true Lie algebra creates difficulties in canonical 



quantum gravity. Isham and Kuchaf |20| found a way in which the Dirac 'algebra' of gravity 
modified by the imposition of Gaussian coordinate conditions can be turned into a true Lie 
algebra. Kuchaf and Torre |H| showed that this procedure is physically equivalent to the 
coupling of gravity to a Gaussian reference fluid that has the structure of a non-rotating, 
heat-conducting dust. We are presently dealing with a simpler and more natural matter 
source, namely, a general incoherent dust; however, the algorithm to be followed is essentially 
the same. The first step is to replace the old constraints (3.20)-(3.21) by an equivalent set of 
constraints based on the resolution of the old constraints with respect to the momenta P and 
Pk conjugate to the time T and frame Z k variables. This task is simplified by the fact that 
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the energy density H±(x) depends on Pk(x) only through the momentum density H^(x). 
Hence, by using the supermomentum constraint (3.21), the super-Hamiltonian constraint 
(3.20) can be cast into the form 

Hi(x) := P(x) + h(x; g ab} p ab ] = , (4.6) 

with 

h(x) = -\fG(x) , (4.7a) 

G(x; g ab ,p ab ] := (Hfix)) 2 - g"*(x)H°(x)H°(x) . (4.7b) 

The choice of the negative square root in Eq. (4.7a) is mandated by the positivity of P{x), 
Eq. (|3.11|) . Inversely, the positivity of P(x) is a consequence of the new constraint (4.6). 
By introducing the inverse Z k {x) to the matrix Z k ta (x), 

Z\ a (x)Z a k (x) = Si , (4.8) 

the supermomentum constraint (3.21) can be rewritten in the form 

H 1k (x) :=H a (x)Z a k (x)=P k (x) + h k (x; T, Z k , g ab ,p ab ] = , (4.9) 
h k (x; T,Z k ,g ab ,p ab ] = Z a k (x)H^(x) + Jc(x)T a (x)Z%(x) , (4.10) 

in which P(x) has been eliminated by the use of Eq. (4.6). 

The old constraints (3.20)-(3.21) and the new constraints (4.6), (4.9) are completely 
equivalent (in the region G(x) > 0, H^(x) < of the gravitational phase space). However, 
the momenta Pk = (P, Pk) in Eqs. (4.6) and (4.9) are separated from the rest of the canonical 
variables. This has an important consequence for the constraint algebra, namely, that the 
Poisson brackets of the new constraint functions H^ K (x) = (H^(x), H^ k (x)) must strongly 
vanish: 



{H u (x),Hi K (x')} = . (4.11) 

This remarkable fact follows by a simple argument [10] that circumvents the tedious algebra: 
Because the new constraints are equivalent to the old constraints and the Poisson brackets 
of the old constraints weakly vanish, the Poisson brackets of the new constraints also must 
weakly vanish. However, because the momenta Pk appear in the constraints without any 
coefficients, these brackets cannot depend on Pk- Thus, the new constraints cannot help 
in any way to turn those brackets into zeros; this implies that the brackets must vanish 
strongly. 

Thus, the new constraints generate an Abelian algebra (4.11). Let us smear these con- 
straints by a set N^ K (x) = (N^(x), N^ k (x)) of externally prescribed scalar functions: 

H[N^] := [ d 3 xN\x)H^(x) , (4.12) 
miV T ] := / d 3 xN ]k {x)H 1k {x) . (4.13) 
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The smeared constraints generate through the Poisson brackets the dynamical changes of 
the field variables T, Z k , P, P k , g ab , and p ab . In particular, H[N^} generates the change 

f(x) := {T(x), H[N 1 }} = N\x) , (4.14a) 
Z k [x) := {Z k (x), H[N^]} = , (4.14b) 

and H[W] generates the change 

f(x) := {T(x), H[N^}} = , (4.15a) 
Z k (x) := {Z k (x), H[N^]} = N^ k (x) , (4.15b) 

of the dust coordinates T(x) and Z k (x). Equations (4.14) tell us that H [iV T ] displaces the 
hypersurface by the proper time N^(x) along the flow lines of the dust (Fig. 2). The 
Poisson bracket of any remaining variable (like g a b{x)) with H[N^] then gives the change 
of that variable under such a displacement. Equation (4.15) tells us that H[W] does not 
change T(x) but generates a displacement to the flow line whose Lagrange coordinates differ 
by the amount N^ k (x) from the Lagrange coordinates z k = Z k (x) of the flow line passing 
through x. This operation does not leave the hypersurface F(E) fixed; rather, it displaces 
each point of the hypersurface along the hypersurfaces of the dust foliation T (Fig. 3). 
Again, the Poisson bracket of any remaining variable with H[N^\ generates the change of 
this variable under such a displacement of the hypersurface. 

The new form of the constraints leads to the familiar form of the Heisenberg equations 
of motion. Let F[Z K , Pk, g,p] be any dynamical variable on the phase space and choose a 
foliation Z^(x) = Z K (t,x) of T x S. The gravitational variables g a b{ x ) an d p ab (x) evolve 
according to the Hamilton equations of motion from given initial data. These equations 
along with the Hamilton equations (4.14) and (4.15) determine g a b(x; Z^} and p ab (x;Z^]; 
the constraints (4.6), (4.9) then give Pr(x; Zf\. In this way the dynamical variable F 
becomes a functional of Zf ', 

F\Zf\ := F [z^P K [Z^g[Z^p[Z^]\ . (4.16) 

The rate of change of this functional with t is given by 

F[Z t K ] = [ £x^p-Z?{x) = [ d 3 x^^-Nj K (x) . (4.17) 



'6Z?(x) ty ' h 5Z t K (x) 

The symbol d/5Z^(x) denotes the total variational derivative of F\Zf\ with respect to 
Z^(x), which takes into account both the explicit dependence of F[Z K , Px,g,p] on Z K and 
the implicit dependence through the dynamical variables g, p, and Pk- On the other hand, 
F is also given by the Poisson bracket of F[Z K , Pk, g, p] with the Hamiltonian H[N^ K \. 
By comparing the two expressions for F for an arbitrary choice of foliation, we obtain the 
many-fingered time Heisenberg equations of motion 

j^ = {F, HM*» = + (4.18) 

Here, h K (x) = (h(x),h k (x)) are the true Hamiltonian and momentum densities from 
Eqs. (4.7) and (4.10). 
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C. Representing LDiff(T x S) 



In Eqs. (4.12)-(4.13), the new constraints were smeared by externally prescribed multi- 
pliers N' (x). The constraints also can be smeared by vector fields V\, V 2 on T x S in 
the following way. Consider V\ and V 2 as elements of LDiff(T x S) under the Lie bracket 
operation — [Vi, V2] which, up to a sign, is given by the commutator. When restricted to 
an embedding Z K (x), these vector fields become dynamical variables on the phase space of 
the dust. Thus, Vi and V 2 can be used to smear the new constraints on S; for example, 

V 1 ^H 1 [V 1 ] := Jj 3 xV?(Z L (x))H 1K (x) . (4.19) 

In this way each vector field on T x S is mapped into a dynamical variable over the phase 
space of dust coupled to gravity. Equations (4.6), (4.9), and (4.11) then guarantee that 

{H^[V 1 ],H^[V 2 ]} = H ] [-[V 1 ,V 2 ]] , (4.20) 

i.e., that the mapping (4.19) is a homomorphism from LDiff(Tx S) into the Poisson bracket 
algebra over the phase space. In this manner, the introduction of dust enables us to represent 
the Lie algebra of the dust spacetime diffeomorphism group by canonical transformations. 

D. New Constraints Lead to a True Algebra for Vacuum Gravity 

The new Hamiltonian constraint (4.6) has a feature that is quite astounding: The Hamil- 
tonian h(x) of Eq. (4.7) does not depend on the dust variables, but only on the geometric 
variables g a b, p ab . We have not encountered such behavior in any other physical system. 
Even for the simplest parametrized systems, like a massless scalar field on a Minkowskian 
background ||23]| , the Hamiltonian depends on the many-fingered Minkowski time. In con- 
trast, the Hamiltonian (4.7) turns out to be independent of the time T(x) and frame Z k (x) 
variables. 

This remarkable feature has an equally remarkable consequence: The Poisson brackets 
{H^(x), H^(x')} = do not have any cross terms between P(x) and h(x), and hence h(x) 
must have vanishing Poisson brackets among themselves: 

{h(x;g,p},h(x';g,p]} = 0. (4.21) 

Since the Hamiltonian h(x) is a square root (4.7), it follows that the expressions G(x) under 
the square root must also have vanishing Poisson brackets: 

{G(x;g,p),G(x';g,p)} = 0. (4.22) 

This property is vital for the definition of the Hamiltonian operator h(x; g,p] by spectral 
analysis (Section 6c). 

We proved Eq. (4.22) by coupling geometry to dust. However, Eq. (4.22) is an identity 
involving only the geometric variables g a b and p ab . Thus, it must hold irrespective of whether 
geometry is coupled to dust, to any other matter system, or whether it is left alone in vacuum. 
Equation (4.22) must simply hold by virtue of the Dirac 'algebra' (4.1)— (4.3) among the 
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gravitational expressions H^(x) and H^(x). We verify this assertion directly in Appendix 
E. 

We have found a simple quadratic combination G(x) of the gravitational super- 
Hamiltonian H±{x) and supermomentum H^(x) that strongly commutes with itself, 
Eq. (4.22). Because G(x) is a scalar density of weight 2 constructed solely from the ge- 
ometric variables, its Poisson bracket with the gravitational supermomentum is given by 

{G(x), H^(x')} = G >a {x)5(x, x) + 2G(x)5 ja (x, x) . (4.23) 

Of course, the gravitational supermomenta represent LDiflE, 

{H% (x), H^(x')} = H^(x)5 ta (x, x') - (ax <-> bx) . (4.24) 

Thus, the geometric expressions G{x) and H^(x) form a true Lie algebra that is a semidirect 
product of the Abelian algebra generated by G(x) and LDiffE generated by H^(x). 

In vacuum gravity, the imposition of the usual super-Hamiltonian and supermomentum 
constraints H^(x) and H^{x) is clearly globally equivalent to the imposition of the new set 
of constraints, G(x) and H^(x): 

Hf{x) = = H°(x) <^> G(x) = = iff (x) . (4.25) 

Therefore, we have succeeded in casting vacuum gravity into a form in which the constraints 
close according to a true Lie algebra. This is something that previously has been achieved 
only for systems in two-dimensional spacetimes, such as strings P4J and 2-dimensional 



gravity |25 



There is only one fly in the ointment. In vacuum gravity, the Hamiltonian vector field 
generated by G(x) vanishes on the constraint surface (4.25). Therefore, G(x) does not 
generate a motion on the constraint surface, and it cannot be used for foliating the constraint 
surface into orbits. Such a difficulty does not exist for gravity coupled to dust: As long as 
there is any dust at x, 

G(x) = P 2 (x) > . (4.26) 

The full significance of the new form of the constraints for vacuum gravity remains to be 
investigated. 



V. CANONICAL DUST SPACE VARIABLES 



A. The Dust Space Variables and the Constraints 

Our description of dust coupled to gravity was based on the conjugate pairs Z k {x) and 
Pk{x), T(x) and P(x), and g a b{x) and p ab (x) of canonical variables. The variables T(x), 
P(x), g a b(x), &iadp ab (x) are (t-dependent) tensor fields on E. The variable Z k [x) is the local 
coordinate chart expression of a (t-dependent) mapping Z : E — > S from the space manifold 
E to the dust space S. Its conjugate Pk(x) is the coordinate expression of a (t-dependent) 
mapping from E to T*S, with density weight 1 on E. The Hamiltonian formulation of 
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the combined system is characterized by the constraints H^(x) = and H^ k (x) = which 
generate the abelian algebra (4.11). 

Canonical quantization is vastly simplified by the introduction of an alternative set of 
canonical variables which reflect the fact that the dust particles naturally define a preferred 
set of spatial coordinates on E. In this section, we construct such variables, and explain 
how they are geared to the abelian algebra (4.11). 

Start from the tensorial variables T(x), P(x), g a b(x), and p ab (x) on E, and let the map- 
ping Z : E — > S induce a corresponding set of tensors T(z), P(z), g ki (z), and p ke (z) on the 
dust space S. In terms of local coordinates, these 5-variables are defined by 



T(z) :- 
P(z) ,- 

p k \z) :-. 



T{X{z)) , 
dX(z 



dz 



P(X(z)) , 



X\ k {z)X\,{z)g ab {X{z)) , 



dX(z) 



dz 



Z k , a (X(z))Z\ b (X(z))p° b (X(z)) 



Here, the t-dependent mapping X is simply the inverse of Z, 

X : S -> E , X := Z- 1 , 



(5.1a) 
(5.1b) 
(5.1c) 
(5.1d) 

(5.2) 



and \dX(z)/dz\ is the Jacobian for the change of variables x a = X a (z). 

The 5-variable T(z) measures proper time along the dust worldline whose Lagrangian 
coordinate label is z k , and P(z) is the dust rest mass on E per unit coordinate cell d 3 z. The 
variables g ke (z) are the components of the metric tensor 



ds 2 = g ke (z) dz k dz l 



(5.3) 



that measures proper distance ds in E between neighboring dust particles with Lagrangian 
coordinate labels z k and z k + dz k . 

We will now show that the iS-variables T(z), P(z), g ke (z), and p ke (z), along with the 
dust frame Z k (x) and the new supermomentum H^ k (x), are canonically conjugate. That is, 
we will argue that the only nonvanishing Poisson brackets among these variables are 



{Z k (x),H^(x')} = 5 k 5(x,x') , 
{T(z),P(z')} = 6(z,z>) , 

P m V)} = \ ( W + 8(z, z>) . 



(5.4a) 
(5.4b) 

(5.4c) 



Recall (Section 4b) that the constraint functions H^ k (x) have strongly vanishing Poisson 
brackets among themselves. They are thus natural candidates for a new set of canoni- 
cal momenta. Equation (5.4a) maintains that their conjugate coordinates are the frame 
variables Z k (x). This is obvious from the construction (4.9) of H^ k (x). 

The Poisson brackets (5.4b) and (5.4c) exemplify a general rule: the mapping E — > S 
turns canonically conjugate fields on E into canonically conjugate fields on S. This is most 
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easily established for a scalar field and its conjugate, such as the proper time (5.1a) and the 
rest mass density (5.1b). In this case, we first write 

T(z) = d 3 x T(x) 5(x - X(z)) , (5.5a) 

P( z ) = d 3 x P{x) 5{Z{x) - z) , (5.5b) 

and then directly verify that the Poisson brackets 



, :; , 5T{z) 5P(z') 5T(z) 5P(z') 



{T W , P(z'» = J s ft {^JJ^f - Jf^^f) + • • ■ < 5 - 6 > 

yield the desired result (5.4b). (The unwritten terms in Eq. (5.6) involve functional deriva- 
tives with respect to Z k (x), Pk{x), g a b(x), and p ab (x), and they all vanish.) It is also clear 
that the Poisson brackets of T(z) among themselves and of P{z) among themselves strongly 
vanish. It is straightforward to extend these arguments to arbitrary canonically conjugate 
tensorial variables and, in particular, to show that g ki (z) and p ke (z) are canonically conju- 
gate. 

The iS-variables T(z), P(z), g M (z), and p kt [z) obviously have vanishing Poisson brack- 
ets with the frame variables Z k (x). The only remaining task is to show that their brackets 
with the new momentum H^k{x) also vanish. To do that, let us recall the relationship 
H]k(x) = Z k (x)H a (x) between the new supermomentum H^k(x) and the original supermo- 
mentum H a (x). Since the matrix Z k depends only on the frame variables Z k (x) and the 
5-variables (5.1) are independent of Pk{x), they have vanishing brackets with Z%(x). Thus, 
the Poisson brackets of the iS-variables with H^{x) are proportional to their brackets with 
H a (x). However, the Poisson brackets of the S- variables with H a (x) must vanish. Recall 
that the Poisson brackets of any dynamical variable with the supermomentum smeared by 
an externally prescribed vector field N a (x) on S, 

H[N] := J d 3 xN a (x)H a (x) , (5.7) 

give the change in that variable under LDiflX. Since the 5-variables are constructed so 
that they are invariant under DiffE, their Poisson brackets with H[N] and hence with H a (x) 
vanish. It follows that their brackets with H^{x) also vanish. This completes our proof that 
the S- variables (5.1) along with the pair Z k (x), H^k(x) form a canonical chart. 

In terms of the new canonical variables T(z), P(z), g k i(z), p ke (z), and Z k (x), H^ k (x), 
the supermomentum constraint is simply the condition that the canonical momentum H^k{x) 
vanishes. The super-Hamiltonian constraint H^(x) = can then be mapped to the dust 
space S, which yields an equivalent constraint 

dX(z) 



H,(z) := 



dz 



H^X(z)) = P{z) + h(z; g ke ,p ke ] = . (5.8) 



Here, the Hamiltonian density h(z) = —^jG(z) is obtained from h(x) of Eq. (4.7) by re- 
placing the E-tensors g a b{x) and p ab (x) with the corresponding iS-tensors gkt{ z ) and p M (z). 
The constraint H f (z) is expressed entirely in terms of the iS-variables. Its Poisson brackets 
with H^ k (x) therefore vanish. We conclude that the constraints H^{z) and H^ k (x) still form 
an abelian algebra. 
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B. Noether Charges 



The Noether charge Q[&] from Eq. (3.25) can be expressed in terms of the new canonical 
variables by a simple change x a = X a (z) of integration variables. This gives 

Q[&] = / d 3 x #(Z(x))P(x) (5.9a) 

= f d 3 z &(z)P(z) . (5.9b) 
J s 

As discussed in Sees. 3b and 2e, Q\&\ generates a canonical transformation that yields 
the change of any dynamical variable induced by the change of the initial hypersurface of 
constant dust time T(y) = t = const. Among the new canonical variables, only the dust 
proper time T(z) is affected by Q[&\: 

{T(z),Q[*]} = *{*) ■ (5-10) 

The canonical transformation generated by Q[&] produces a translation of the proper time 
r along the dust worldline z k by the amount "d(z). 

— * 

The expressions (4.6)-(4.10) can be used to write the Noether charge Q[£] from Eq. 
(3.24) in terms of the new canonical variables. The result is 

Q[£\ = - J^d 3 x e(z(x))P k (x) (5.11a) 
= - jj 3 xe(z{x))H^ k {x) + J^z £ k (z)(H°(z) + P(z)T k (zf) , (5.11b) 

where H k (z) is obtained by mapping the gravitational supermomentum H^{x) from space 
E to the dust space S: 



H G k {z) 



dX(z) 



dz 



X\ k {z)H<;{X{z)) . (5.12) 



Thus, H)*(z) depends only on g ki (z) and p ke (z). 

The Noether charge Q[£] can be recognized as the canonical generator of dust space 

— * 

diffeomorphisms, DifLS. That is, for any canonical variable F, the Poisson brackets {F, Q[£\} 
give the change in F due to an infinitesimal diffeomorphism of S generated by the vector 
field £(z). If F is one of the 5-tensors T, P, g, or p, then the Poisson brackets {F, Q[£\} 
equal the Lie derivative £^F. For the dust frame variable Z k [x) and its conjugate H^ k (x), 

the transformation generated by Q[£] is 

{Z k (x),Qi£]} = -e(Z(x)) , (5.13a) 
{Hi k (x),Q$} = ^,k(Z(x))H n (x) . (5.13b) 

The charges Q[&] and Q[£\ close under the Poisson brackets according to the relations (3.28). 

As noted in Sec. 3b, the Noether charges Q{&\ and Q[q\ have vanishing Poisson brackets 
with the constraints H±(x) and H a (x). Since the new set of constraints H^ k (x) = and 
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H^(z) = are equivalent to the old set, the Poisson brackets of the charges with the new 
constraints should vanish at least weakly, that is, when the constraints hold. Indeed, the 
Poisson brackets of Q[&] with H^ k (x) and H-\(z) vanish strongly. The Poisson brackets of 
Q[$] with H^ k {x) vanish when the constraints hold, as seen from Eq. (5.13b). The remaining 
Poisson brackets are 

{H 1 (z),Q[Z]}=(H 1 (z)e(z))^ , (5.14) 

and they also weakly vanish. Expression (5.14) shows that H^-(z) transforms as a scalar 
density under DifLS. 



VI. QUANTIZATION 
A. Operator Constraints and the Functional Schrodinger Equation 

In the Dirac method of quantization, constraints are turned into operators and imposed 
as restrictions on the state functionals of the system. With the standard form (3.20) of 
the super-Hamiltonian constraints, this procedure leads to the Wheeler-DeWitt equation 
for gravity coupled to a source. This is a second-order variational differential equation in 
the metric variables whose solution space carries no obvious Hilbert space structure. If the 
source happens to be dust, there is an additional difficulty: the square-root Hamiltonian 
density (3.19) is difficult to define by spectral analysis because the expressions under the 
square root do not commute, Eq. (4.5). 

In order for dust to serve as a clock, the super-Hamiltonian constraint should be resolved 
with respect to the momentum conjugate to the dust time variable prior to quantization. 
The constraint H^(x), or equivalently H^z), has this desired form. We choose to work 
primarily with the dust-space variables, so the quantum states of the system are functionals 
&[Z,T,g] of the canonical coordinates Z k (x), T(z), and g ki (z). From the expression (5.8) 
of the classical constraint H^(z) = 0, the operator condition 

H^z;P,g,p]^[Z,T,g] = (6.1) 

yields the functional Schrodinger equation 

.mZ,T,g] 



5T(z) 



= h(z;g,p}*[Z,T,g} . (6.2) 



In addition to Eq. (6.2), the state functionals must obey the operator version of the classical 
constraint H^ k (x) = 0, namely, 

H ]k {x)V[Z,T,g\ = . (6.3) 
Because H^ k {x) is canonically conjugate to Z k (x), H^ k {x) is the momentum operator 
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Equation (6.3) then implies that the quantum states must be independent of the dust frame 
Z k {x). Therefore the state functionals *f?[T, g] satisfy the functional Schrodinger equation 

5^\T a] 

i -sTW = h ^^p]^[T,g], (6-5) 

with no further restrictions. 

The functional differential equations (6.1), (6.3) are consistent only if the commutators 
of the operators H^{z) and H^ k (x) replicate the classical Poisson bracket algebra 

{H r (z),H ] (z')} = 0, (6.6a) 
{Hi(z),H Vc (x)} = , (6.6b) 
{H^(x),H n (x')} = . (6.6c) 

With the representation (6.4) for the momentum operator the subalgebra (6.6c) is indeed 
carried over into the quantum theory. The brackets (6.6b) are also satisfied quantum me- 
chanically since each of the operators P, g, p that appear in H^{z) commutes with H^{x). 
We shall proceed under the assumption that there exists a factor ordering and regulariza- 
tion of h(z; g,p] such that the quantum operators H^{z) commute among themselves, and 
thereby replicate the vanishing of the classical Poisson brackets (6.6a). 



B. Time Independent Schrodinger Equation 

H-\{z) does not depend on T(z) so that 

[H,(z),P(z>)} = 0. (6.7) 

Furthermore, the operators P(z) commute among themselves, and the super-Hamiltonian 
operators H^(z) commute by virtue of Eq. (6.6a). It follows that the Schrodinger equation 
(6.5) and the eigenvalue equation 

P{z) * = P{z) * (6.8) 

have common solutions ^[P; T, g\. These are functionals of T(z) and g^eiz) labeled by the 
eigenvalues P(z) of the momentum operator P(z). The general solution of Eq. (6.8) has 
the form 

*[P; T,g] =iP[g]exp(iJ s d 3 zP(z)T(zf) . (6.9) 

By substituting expression (6.9) into Eq. (6.5) we learn that i/)[g] must satisfy the time- 
independent functional Schrodinger equation 

h(z;g,p]^[g] = -P(z)^[g] . (6.10) 

Equations (6.9) and (6.10) accomplish the separation of the many-fingered time T(z) from 
the metric g ke (z). 
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C. Spectral Analysis 



The classical Hamiltonian density h(z) is the dust space counterpart of the square root 
Hamiltonian density h(x) = —^jG(x), as defined in Eq. (4.7). The square root structure of 
h(x) is inherited by h(z), so that h(z) = —^jG(z) where 

G(z; g,p] := (H G ± (z)f - g k \z)H G (z)Hf(z) . (6.11) 

As a consequence, the operators h(z) that appear in the functional Schrodinger equation 
must be defined by spectral analysis. The first step is to turn the classical expressions (6.11) 
into operators G(z; g,p]. The assumption that the operators H-\(z) commute translates 
into the assumption that there exists a factor ordering of g and p, and a regularization of 
G{z), such that 

[G(z;g,p],G(z';g,p]] = 0. (6.12) 
The operators G(z) then have common eigenfunctionals <j>[G; g], 

G(z) cf>[G; g] = G(z) <f>[G; g] , (6.13) 

with eigenvalues G(z). 

The classical variables (6.11) are not positive definite, so the eigenvalues G(z) are not 
necessarily positive everywhere on S. Let us define a Hilbert space 7i + that is spanned by 
the eigenfunctionals of G(z) with positive eigenvalues G(z) > 0. On H + , we can define the 
operator 

h{z) = -\[G~{z) (6.14) 
by the requirement that it have eigenfunctionals <p[G; g] with eigenvalues 

h{z) = -^G{z). (6.15) 

Because h(z) is the Hamiltonian of the system, a state functional in 7i + will remain in 7i + 
throughout its dynamical evolution. 

If an eigenfunctional <j)\G\ g] with eigenvalue G(z) is known, then infinitely many eigen- 
functionals corresponding to other eigenvalues can be generated by the action of dust space 
diffeomorphisms S e Diff«S. This result is expressed by the following theorem, which is 
proved in Appendix G: Let <f>[G; g] be an eigenfunctional of G(z) with eigenvalue G(z). 
Then <f)[G; H*gr] is an eigenfunctional of G(z) with eigenvalue (S _1 )*G(z). 

Disregarding global difficulties, a given dust space scalar density G(z) > of weight 
2 can be transformed by the action of an appropriate diffeomorphism S _1 into any other 
arbitrary function G*(z) > 0. In particular, it can be transformed into G*(z) = 1. This 
means that G(z) has a continuous spectrum, i.e., that there is an eigenfunctional <p[G; g] 
for any G(z) > 0. 
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These arguments also indicate that the spectrum of the collection z G S of commuting 
operators G(z) is infinitely degenerate. Indeed, let S G Diff iS be an arbitrary volume- 
preserving diffeomorphism, det(d'H.Q(z)/dz) = 1. For G(z) = 1, we have also (Hq 1 )*G(^) = 
1. Therefore, according to the theorem stated above, if 4>[l;g] is an eigenfunctional of 
G(z) with eigenvalue 1, then 0[1; H^g] is another eigenfunctional of G(z) with eigenvalue 
1. This shows that there are at least as many eigenfunctionals of G(z) with eigenvalue 1 as 
there are volume-preserving diffeomorphisms. When we act on the argument of this set of 
eigenfunctionals by an arbitrary diffeomorphism H, we turn it into a set of eigenfunctionals 

S*S* g] , H G Diff <S (6.16) 

where each eigenfunctional has the eigenvalue G(z) := (det(dS(z) / dz))~ 2 . In this way, we 
conclude that an arbitrary level G(z) > is infinitely degenerate. 

It remains an open question whether there is more degeneracy of the G(z) > levels than 
that corresponding to the action of volume-preserving diffeomorphisms, and, if so, how the 
additional degeneracy should be labeled. We do not need to answer this question. Formally, 
the degeneracy of the levels G(z) > can be taken into account by the introduction of 
a degeneracy index d. The Hilbert space H + is spanned by the eigenfunctionals <f> d [G;g] 
which satisfy 

G(z) 4> d [G; g] = G(z) 4> d [G; g] ■ (6.17) 

An arbitrary state ^[g] in 7i + can be written as a linear superposition 

*bl=EL DG C d [G] 4> d [G] g] (6.18) 

d JG{z)>Q 

of the eigenfunctionals </> d [G; g] with coefficients C d [G\. These eigenfunctionals form a basis 
for the energy-density representation (6.14)-(6.15). 

From the considerations above, the general solution of the functional Schrodinger equa- 
tion (6.5) can be written in the form 

^[T,g] = U[T-T ; g,p]*[T ,g} , (6.19) 

where 

U[T-T ; g,p] := exp(-i J s #z{T(z) - T (z))h(zf) . (6.20) 

Here, *S?[T ,g] is an arbitrary initial state in H + , and the operator h(z) is defined by 
spectral analysis. The state Vl/[T, g] at the dust time T(z) follows from the initial state 
Vl/[T , g] at the dust time T (z) by the action of the many-fingered-time evolution operator 
U[T-T ;g,p}. 



D. The Action of Symmetries on States 

The wave functional *f?[Z, T,g] characterizes the quantum state of the gravitational field 
on the hypersurface specified by the dust time T(z). Such a description is not unique. First 
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of all, the dust time is not unique because the time foliation T depends on the choice of the 
fiducial hypersurface T(z) = 0. By changing this hypersurface, the time variable T(z) is 
changed according to 

T'(z) = T(z) + &(z) . (6.21) 

(This is Eq. (2.27) expressed in the dust coordinates.) The variables T(z) and T'(z) con- 
nected by Eq. (6.21) specify a given hypersurface relative to two choices of dust time, T and 
T . Thus, the fiducial hypersurface T(z) = of the time foliation T is also the hypersurface 
T'(z) = @(z) of the time foliation T. 

The wave functional &[Z, T, g] does not provide a unique description of the gravitational 
field on a given hypersurface for a second reason: The Lagrangian coordinates z k that label 
the dust worldlines are not unique. A relabeling of the dust worldlines is brought about by 
a diffeomorphism 3 : S — > S of the dust space, as discussed in Sees. 2e and 5b. Under the 
action of S, the dust time T and the metric g change by the pullback mappings 

T':=ToS, g':=E*g. (6.22) 

According to Eq. (2.26), the action of a diffeomorphism S on the dust frame Z is given by 
Z'-.^S^oZ. 

The value of the state functional should not depend on the specification of the hyper- 
surface. Thus, if *&[Z, T, g] describes the state of the gravitational field relative to the dust 
time T and *&'[Z, T' , g] describes the same state relative to the dust time T', then these 
two wave functionals should be equal: 

T' = T + 0, g] = y[Z,T,g] . (6.23) 

Likewise, the value of the state functional should not depend on the choice of Lagrangian 
coordinate labels for the dust particles. Therefore, 

V\Z' = '3r x oZ, T' = ToS, g' = E*g] = ^[X,T,g] (6.24) 

where *&[Z, T, g] describes the state of the gravitational field with respect to the Lagrangian 
coordinates z k and Sf?'[Z',T',g'] describes the same state with respect to the Lagrangian 
coordinates z k . 

These changes in the description of the state are generated by the Noether charges. Thus, 
the state Vl/' from Eq. (6.23) can be obtained from the state ^ by the action of the operator 

U[S] :=exp(-iQ[0]) (6.25) 
generated by the Noether charge 

Q[0] = f d 3 z &{z)P(z) . (6.26) 

Indeed, we have 

T', g] = *[Z, T' - 0, g] = U[@] *[Z, T', g] . (6.27) 
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Similarly, let 3 fc (-2, cr) denote a one-parameter group of diffeomorphisms with E k (z, 0) = z k 
corresponding to the identify and S k (z, 1) = S k (z) corresponding to a given final diffeomor- 
phism. As usual, define the infinitesimal generators 



-/■,.. . d~ k (z,a) 



£(z) ■= ^ ' (6-28) 

and the Noether charges Q[£ a ] corresponding to the classical charges (5.11). Then the state 
from Eq. (6.24) can be obtained from the state 

*'[Z',T',g'] = *[~oZ', T'oE" 1 , (H" 1 ) *g'\ = U[~] *[Z>, T', g'\ , (6.29) 

by the action of the operator 

U[S] : = T CT exp(-* £ da • (6-30) 

Here, T CT demands that the operators Q[£ a ] be ordered so that Q[CrJ precedes ^ 
o"2 > cr 1 . 

The consistency of the formalism requires that physical states be transformed by the 
change of description again into physical states. Thus, if ^ satisfies the constraints (6.1) 
and (6.3), the state Vl/' defined by Eq. (6.23) should also satisfy the constraints. This follows 
from the commutator of the Noether charge (6.26) with the constraint operators: 

[Q[0], ff t (z)] = = [Q[&], H lk (x)\ . (6.31) 

Equations (6.27) and (6.31) imply 

H^z) *'[Z,T,g] = U[e]H,(z) *[Z,T,g] , (6.32) 

and similarly for H^ k (x). Hence, if ^ satisfies the constraints, VP' also satisfies the con- 
straints. 

Similarly, if \I> satisfies the constraints, then ty' as defined in Eq. (6.24) should also 
satisfy them. This follows from the commutation relations 

[Hl k (x),Q[£]}=ie, k (Z(x))H n (x) , (6.33a) 
[Hi(z),Q[£]]=i(Hi(z)e(z)), k , (6.33b) 

which are the quantum counterparts of the classical Poisson brackets (5.13b), (5.14) between 

— # 

Q[£\ and the constraints H^ k (x), H^{z). Equations (6.33) show that if *f?[Z, T, g] satisfies the 
constraints, then Q[£\ &[Z,T,g] also satisfies the constraints. In turn, if Q[£] *&[Z,T,g] 
satisfies the constraints, then *S?'[Z,T,g] = U[S]^f[Z,T, g] also satisfies the constraints. 
We conclude that the description of a physical state can be changed by the action of the 
operators U[S] and U[S] generated by the Noether charges. 
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E. The Inner Product 

The transition from the Hamiltonian constraint H±(x) to the new constraint H-^(z) leads, 
upon quantization, to the replacement of the Wheeler-DeWitt equation by the functional 
Schrodinger equation (6.5). Unlike the Wheeler-DeWitt equation, the Schrodinger equation 
has (at least formally) a conserved positive-definite inner product || . This, of course, is the 
main motivation behind rewriting the Hamiltonian constraint in the forms (4.6) and (5.8). 

Let *f?[g] be a functional of Riemannian metrics g(z) G RienuS in the dust space S. The 
inner product between two such functionals is defined by the functional integral 

(*i|* 2 > := / Bg*i\g]*2\g] ■ (6-34) 

J RiemS 

Those functionals *f?[g] that are square integrable under the product (6.34) form a Hilbert 
space TC. 

The functionals &[g] can be interpreted as states of the gravitational field on a fixed 
initial hypersurface T (z): 

*[T ,g] :=*[<?] . (6.35) 

By virtue of their independence from Z, the state functionals *&[Tq, g] automatically satisfy 
the operator constraints (6.3). Let us evolve these states by the Schrodinger equation (6.5). 
The solutions Vf^T, g], from Eq. (6.19), form a space which we call Tip By fixing T(z) to be 
the given initial hypersurface T (z), a functional *&[T, g] from is mapped to a functional 
(6.35) from H. 

We now ask whether the inner product 

(*i|* 2 ):=/ Dg* 1 [T,g]* 2 [T,g] (6.36) 

J RiemS 

between two states from the solution space TC^ depends on the choice of hypersurface T(z). 
Formally, the answer is no, provided that the Hamiltonian density h(z\ g,p] of Eq. (6.5) 
is a self-adjoint operator under the inner product (6.34). However, h is not self-adjoint 
on the whole space 7i, but only on the linear subspace 7i + C Ti which is spanned by the 
simultaneous eigenfunctions (6.13) of the operators G(z) with positive eigenvalues G(z) > 0. 
Hence, if the inner product (6.36) is to be independent of the choice of hypersurface, we 
must limit the solution space Hi to Tit . 

To summarize, the state of the gravitational field can be represented either as an element 
of the space of solutions or as an element of the space of initial states 7i + . Conservation 
of the inner product guarantees that the value obtained for the inner product is the same 
irrespective of the representation. 

F. Observables 

Any functional F[g,p] of the metric g(z) and extrinsic curvature p(z) represents a grav- 
itational dynamical variable that is invariant under DiffS. We will call the space of such 
variables T . If the dust particles are labeled by a fixed set of Lagrangian coordinates z k , 
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then in principle the full metric g^{z) on a given hypersurface, not just the intrinsic ge- 
ometry of the hypersurface, can be measured. A corresponding statement holds for the 
extrinsic curvature. If we want, we can restrict our attention to variables in T that do not 
depend on the choice of Lagrangian coordinates for the dust; such variables have vanishing 
Poisson brackets with the generator Q[£] of dust space diffeomorphisms. In the following, it 
makes no difference whether we take the position that all the dynamical variables in T are 
measurable, or only those that are invariant under DifLS. Thus, we leave open the issue of 
which class of variables should be called observables. 

In quantum theory, an observable from T is represented by an operator 

F = F[g,p] (6.37) 

that acts on 7i and is self-adjoint under the inner product (6.34). The statistical interpre- 
tation of quantum geometrodynamics is based on two fundamental assumptions. The first 
one is the mean value formula: When F is measured on an ensemble of systems described 
by the state functional &[T ,g] G H on a hypersurface T (z), the mean value of the result 
is given by the expression 

(F) = , (6.38) 

with the inner product (6.34). By the standard argument, it follows that the permissible 
values of F are the eigenvalues of F and, if the system is in one of the eigenstates of 
F, the measurement of F on the hypersurface T (z) yields the corresponding eigenvalue 
with certainty. The second assumption is the reduction postulate: If F is measured on a 
hypersurface Tq(z) and those systems for which the measurement yields the eigenvalue F' 
are collected into a new ensemble, then the state functional *&'[T ,g] that describes the 
new ensemble is the orthogonal projection (under the inner product (6.34)) of the original 
state *[T , g] to the linear subspace of H spanned by the eigenstates of F belonging to the 
eigenvalue F' . 

So far, observables have been represented by operators acting on a state ty[To,g] e T~C 
at the time of measurement T (z). The functional 

<f>[T ,g] = F*[T ,g] (6.39) 

lies in the same instantaneous Hilbert space Ti as the original functional SI/. Let us consider 
the functionals \I/ and <fr as initial data and evolve them by the Schrodinger equation (6.5) 
away from the initial hypersurface. In this way, we obtain states (6.19) that belong to the 
solution space Tlf. 

*lT,g] = U[T-T ; g,p] *[T ,g] , (6.40) 
*[T,g]=U[T-T 0] g,p] <S>[T ,g] . (6.41) 

The state (6.41) is obtained from the state (6.40) by the action of the operator 

F T [T; g,p] = U[T-T ; g,p] F[g,p] U'^T - T ; g,p] . (6.42) 

Equation (6.42) represents the observable F[g,p] which is measured on the hypersurface 
T (z) as an operator Ft on the solution space 7i^. The operator F To is constructed so 
that it commutes with the Hamiltonian constraint H^(z) : 
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[F To ,H,{z)]=0 



(6.43) 



This property guarantees that the action of F To on a state in again yields a state in H^. 

We must now translate into the quantum theory the positivity condition G(z; g,p] > 
on the domain of classical observables F[g,p]. Quantum mechanically, the action of an 
operator F[g, p] on a state *f?[g] G 7i + from the subspace H + spanned by the eigenfunctionals 
of the operators G(z; g,p] with positive eigenvalues G(z) > should leave the state in 7i + . 
If an observable F G T does not have this property, the reduction postulate leads us into 
difficulty The reduced state *&'[T ,g] obtained after the measurement of such an observable 
would not lie in H + . Therefore, if we measure G(z) on the reduced state there would be 
a nonvanishing probability of finding G(z) in the unphysical region G(z) < 0. Furthermore, 
by allowing the reduced state *&'[T ,g] to propagate to a later time T(z) > T (z) via the 
evolution operator U[T— T ; g, p], we would find that the norm (*f?'[T, g]\^'[T, g]) depends 
on T(z) and the total probability is not conserved. 

These considerations force us to define observables as self-adjoint operators F G JF+ on 
7i + rather than on 7i. That is, an operator F is considered to be an observable only if its 
action on the states St from 7i + again yield states in 7i + . Correspondingly, the observables 
F T[) are self-adjoint operators on the space of solutions Tit. Thus, for St an element of Tit, 
the state F To *& also is an element of . This implies that the mean value of F To , defined 
by 

(F To ) = (*|F To *) , (6.44) 

can be evaluated on any hypersurface, not only at the instant T (z) of the measurement. 

The condition for an operator defined on Tt to be an observable can be expressed in 
terms of the projection operator 

V+ : = Yl 0(G(z)) (6.45) 

from TL to TL + . The projector V + is a product of step functions 9 of the commuting operators 
G(z). The Hilbert space H + is spanned by states of the form V + ^, where * G H. An 
operator F on H is an observable if it satisfies the condition 

[F,V + ]V + = (1 -V + )FV + = . (6.46) 

Thus, any variable F[g,p] that commutes with G(z), i.e., that is a constant of motion, is 
automatically an observable. However, F does not need to be a constant of motion to satisfy 
Eq. (6.46); observables do not need to be constants of motion. 

The fundamental problem in interpreting canonical quantum gravity through its coupling 
to dust is the construction of meaningful gravitational observables. Ultimately, one would 
like to exhibit complete sets of commuting observables. 

On our way, we met some gravitational observables. The operators G(z; g, p], equivalent 
to the energy densities h(z; g, p], form a commuting set of gravitational observables on 
7i + . The Noether charges Q[£] and Q[@] represent observables on 7t^ . By using the 
constraints, we can eliminate the embedding momenta and express the charges in terms of 
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the gravitational variables. There are thus some observables we can measure and subject to 
statistical analysis. All of these are constants of motion. 

It is disappointing that the fundamental gravitational variables g(z) = g(z) x and p(z) = 
—i5/5g(z) are not gravitational observables. Although the condition (6.46) is cumbersome to 
handle explicitly, it is nevertheless fairly obvious that g(z) and p{z) do not necessarily leave 
the states in 7i + . The underlying difficulty is that the gravitational Hamiltonian h(x; g,p] 
is not real by its construction on the entire Hilbert space 7i which carries the representation 
of the fundamental operators g(z) and p(z). 

One can illustrate these problems on simple finite-dimensional models, like a relativistic 
particle in the proper-time formalism (Appendix G). In particular, one can explicitly verify 
that even on a flat Minkowskian background the multiplication position operator is not an 
observable. One may hope to clarify the status of fundamental variables on such simple 
models. For the time being, their relation to quantum observables remains unclear. 

G. Quantum Theory in Terms of State Functionals ty[Z,T,g] 

The quantum theory of gravity coupled to dust, as developed thus far, can be expressed 
in terms of the original set of canonical variables Z k (x), Pk(x), T(x), P(x), g a b(x), an d 
p ab (x). In this case, the quantum states are functionals ^f[Z,T,g] that satisfy the operator 
constraint equations 

Hi(x; P,g ab ,p ab ]nZ,T,g] = 0, (6.47a) 
H 1k (x; Z k ,P k ,f,g ab ,p ab ] *[Z,T,g] = . (6.47b) 

Let us assume that Vl/[T, g] is a solution to the functional Schrodinger equation (6.1). Recall 
that such a functional automatically satisfies the operator constraint (6.3) by virtue of its 
lack of dependence on Z k {x). Now choose a mapping Z : £ — > S that associates the points 
in space with particles of dust, and consider the functional defined by 

*[Z,T,g] := ^[(Z-'YT, {Z^fg] . (6.48) 

Here, T = {Z^ V )*T and g = {Z^ v )*g denote the pullbacks by X = Z^ 1 of the space 
tensors T and g to the dust space S. As we show below, the functional (6.48) satisfies the 
constraints (6.47) and represents the quantum state of the system in terms of the original 
canonical coordinates. 

Consider first the supermomentum operator constraint (6.47b). Since H^{x) = 
Z k (x)H a (x), this equation implies H a (x) ^ = which in turn implies that ^ is invari- 
ant under spatial diffeomorphisms. Indeed, let ip c : £ — > £ denote a one-parameter family 
of diffeomorphisms generated by the vector field N(x). Assuming an appropriate factor 
ordering, we have 

= i [ d 3 xN a (x)H a (x) V[Z,T,g] = (6.49) 

cr=0 JT, 

for any wave functional that satisfies the supermomentum operator constraint. By exponen- 
tiation in a, this result yields 



—y[cplZ,cp* a T,cp* a g} 
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*[ ( p*Z, ( p*T,<p*g]=*[Z,T,g] , (6.50) 

and shows that \I> is invariant under finite diffeomorphisms of S. Conversely, a diffeomor- 
phism invariant functional (6.50) will satisfy the supermomentum constraint (6.47b). The 
functional \1/ defined in Eq. (6.48) is just such an invariant functional. This follows from the 
fact that (Z~ X )*T and (Z~ l )*g are invariant under the action of DiffS: For (Z~ V )*T, we 
have 

((</?* Zy 1 )* \ip*T) = (T o (p) o (Z o ip)- 1 = ToZ- 1 

= (Z' l )*T , (6.51) 

and similarly for (Z~ r )*g. Therefore the state functional ty[Z,T,g] defined in Eq. (6.48) 
satisfies the supermomentum constraint (6.47b). 

The functional (6.48), which we assume satisfies the functional Schrodinger equation 
(6.1), also satisfies the functional Schrodinger equation (6.47a). In order to understand this 
result, it is only necessary to clarify the relationship between the operators P, g, and p and 
their dust space counterparts P, g, and p. We will assume that these sets of operators are 
related by the familiar mapping of tensor fields by Z. For example, the action of P(x) and 
the action of P(z) are related by 



P(x) *[(Z- l YT, (Z- 1 Yg]= (P( z )Hf[T,g] 



. (6.52) 

--Z(x), T=(Z-i)T, g=(^- 1 )*s 



There are no factor ordering ambiguities in this relationship. Using an abbreviated notation, 
we can write P(x) = \dZ(x) / dx\P(Z(x)). Observe that if the operator P(x) is represented 
by —i5/5T(x) then P(z) is represented by —iS/ST(z). A relationship analogous to that 
expressed in Eq. (6.52) can be used to define a correspondence between functions of the 
operators P, g, and p, and functions of the operators P, g, and p. In particular, the new 
super-Hamiltonian operator is (using the abbreviated notation) 



H-\(x; P,g,p\ 



dZ(x) 



Hi(Z(x); P,g,p] . (6.53) 



dx 

From this relationship we obtain the correspondence 

Hi(x)*[(Z- 1 yT,(Z- 1 )*g] = ^ H^z)^[T,g] = . (6.54) 

Therefore, if &[T,g] satisfies the Schrodinger equation (6.1), then ty[Z,T,g] as defined in 
Eq. (6.48) satisfies the Schrodinger equation (6.47a). 

The preceding arguments show that the state functional ^[T, g\ can be expressed as a 
functional \& of T, g, and a mapping Z : S — > S through definition (6.48). Different choices 
of mappings Z lead to different functionals however, these functionals are related by 
the unitary transformation U[S] from Eq. (6.30) and represent the same physical quantum 
state. 

The spectral analysis, the construction of the Hilbert space and inner product, and 
the definition of observables can be translated into the language of the original canonical 
variables in a straightforward manner. In particular, the Hilbert space H + is spanned by 
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the eigenfunctionals <j>\{Z 1 )*G; [Z x Yg\ of G(x) with positive eigenvalues G(x) > 0. The 
inner product between two states of the form (6.48) is 

(*!|* 2 > := / Dg tip-^r, g] ^[(Z^fT, g] (6.55a) 

J RiemS 

= / Dg^ 1 [Z,T,g}^ 2 [Z,T,g} . (6.55b) 

Finally, the gravitational observables are operators F[Z,g,p] that correspond, in the sense 
of Eq. (6.52), to operators F[g,p]. Their action on a state vector in the Hilbert space 7i + 
again yields a state in TL + . 



VII. NON ROTATING DUST 
A. Rotation and Vorticity 

In general, no dust time T is better than any other. However, if the dust does not 
rotate as it moves, so that the four-velocity field is hypersurface orthogonal, then there is 
a privileged time foliation 71, namely, the foliation by hypersurfaces orthogonal to the flow 
lines. For such a foliation, Wk = 0. 

The rotation of the dust is measured by the vorticity 2-form 

:= (<*2 + WU a ){8l + U s Up) U b>s] . (7.1) 

For Q a/ 3 = 0, the vector field U a is hypersurface orthogonal. The projectors in Eq. (7.1) 
ensure that flapU^ = 0; hence, VL a p is completely characterized by its orthogonal projection 
Qke '■= ^a / gT a ) / c T /3 ! £. The equations of motion (2.17)-(2.18) imply that 

ft ke = W [kA . (7.2) 

Like W k itself, is conserved along the flow lines of U a : 

£ u n M = n M>a u a = o . (7.3) 

Hence, if the dust is irrotational at a spacetime point y, then the equations of motion ensure 
that it is irrotational for all points along the flow line that passes through y. 

The projected vorticity (7.2) is an exact two-form in the dust space S. Moreover, it 
is invariant under the symmetry transformation (2.27). If Qki vanishes, the velocity field 
Wk is closed; thus Wk can be (at least locally) transformed to zero by the symmetry 
transformation (2.27). This transformation brings us from a foliation T to the privileged 
time foliation 71. 

The foliation 71 also can be constructed by choosing an initial hypersurface orthogonal 
to the flow lines and then Lie propagating it along U a . Because U a is a geodesic vector 
field, Eq. (2.25), the propagated hypersurface stays orthogonal to U a . This amounts to the 
familiar construction of a Gaussian coordinate system [p2 |. The coordinate system (r,z k ) 
based on the time 71 and space S manifolds of an irrotational dust is a Gaussian coordinate 
system. 
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B. Canonical Description of Non Rotating Dust 

The pullback of the vorticity 2-form (7.2) from S to E is 

Qab '■= ^lkeZ k A Z £ >6 = (W k Z k j[0 ) )6] 

= ^ • (7-4) 

Note that in general the direct projection of the spacetime vorticity tensor (7.1) onto £ 
differs from the spatial tensor Vt ab by a term A [a C4], where A a is the spatial projection 
of the dust acceleration A a = U^'VpUa. When the Euler equations of motion hold, the 
acceleration (2.25) vanishes and the spatial projection of the spacetime vorticity coincides 
with the vorticity Q ab of Eq. (7.4). 

Under spatial diffeomorphisms, fl ab (x) transforms as a tensor; this determines its Pois- 
son bracket with the supermomentum H®(x). The Poisson bracket of fl a b(x) and Q c d(x') 
vanishes: 



{n ab (x),n cd (x')} = o . (7.5) 

Finally, the Poisson bracket of fl a b{x) and the super-Hamiltonian is 

ttDc( t\ 

{Q ab (x),H°(x')} = -2—^J-Q c[a (x')S, b] (x,x') . (7.6) 

n ± [x ) 

This result is obtained as follows. Because neither Z k [x) nor W k {x) = —P k (x)/P(x) depend 
on T(x), they have vanishing Poisson brackets with P(x'). The only part of the super- 
Hamiltonian H^(x') that affects these variables is the dust momentum H^(x') in the square- 
root energy-density (3.19). The momentum H®(x f ) generates the transformation of Z k (x) 
and Wk(x) as scalars under DiffE. Therefore, 

{W k (x),H°(x')} = W k , b (x)S b (x)S(x,x') , (7.7a) 

{Z k (x),H° (x')} = Z k , b (x)S b (x)8(x,x') , (7.7b) 

with 

- • (7.8) 

Equations (7.7) allow us to calculate the Poisson bracket 

{W k (x)Z k , a (x),H°(x')} = 2Q ab (x')S b (x')S(x,x') 

+ (u k (x)Z k >b (x)S b (x)5(x,x')) a , (7.9) 

where Vl ab is given by Eq. (7.4). By taking the exterior derivative of Eq. (7.9), we obtain 
Eq. (7.6) above. 

The results (7.5)-(7.6) show that if a new constraint 

n ab (x) = (7.10) 
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is adjoined to the super-Hamiltonian and supermomentum constraints (3.20) and (3.21), 
then the enlarged system of constraints is first class. Moreover, from Eq. (3.14) we have 

Ua = (7.11) 

It is easy to see that both P and are invariant under the symmetry transformations, 
i.e., that their Poisson brackets with the generators (3.24) and (3.25) vanish. Therefore U a 
and fl a b are also invariant under the symmetries: 

{Q ab (x),Qm = = {Q ab (x),Q[i]} . (7.12) 

Thus, the new constraint (7.10) is invariant under the symmetry transformations (3.24), 
(3.25). The condition (7.10) constrains the dust into irrotational flow. The first-class 
nature of the enlarged system of constraints ensures that if the dust is irrotational on an 
initial hypersurface, its dynamics prevents the evolution of vorticity. 

Let us study the vorticity potential Wk(x) := —P k (x)/P(x) from the canonical point of 
view. The Poisson bracket of W k (x) with the super-Hamiltonian H±(x) is different from 
zero; on the other hand, its Poisson bracket with the new constraint H^{x) vanishes: 

{W k (x),H^x')} = . (7.13) 

This reflects the fact that Wk(x) does not change along the worldlines of the dust although it 
does change along the normal direction to an embedding. We can now solve the Heisenberg 
equations of motion (4.18) for the dynamical variables Wk(x) . Any such solution is a 
functional W k (x; Z L ] of the embedding Z L = (T(x), Z e (x)) , as expressed in Eq. (4.16). 
When we substitute Eq. (7.13) into the T-component of the Heisenberg equations of motion, 
we learn that W k (Z L (x)) cannot depend on T(x), i.e., that W k = W k (x; Z ] . Since W k (x) := 
—Pk(x)/P(x) is a spatial scalar, it follows that 

{W k (x),H a (x')} = W k>a (x)S(x,x') . (7.14) 

The definition (4.9) of H^ k {x) then implies that weakly, modulo the supermomentum con- 
straint, 

{W k (x),H ie (x')} = W k , a {x)Zt{x)5{x,x') . (7.15) 
The ^-component of the Heisenberg equations of motion (4.18) thus reads 

dW k (x; Z m 



5Z e (x>) 



Wk,a{x)Z${x)8{x,x') . (7.16) 



The delta function on the right-hand side of Eq. (7.16) tells us that W k (x; Z £ ] is an ultralocal 
functional of Z e (x) , i.e., that 

W k (x; Z L ] = W k (Z\x)) . (7.17) 

It is easy to see that the expression (7.17) solves the Heisenberg equations of motion (7.16). 
In this roundabout way, we recover the canonical counterpart of the spacetime equation 
(2.18). 
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Now impose the condition that the dust does not rotate, so that its vorticity (7.4) 
vanishes. This implies that (locally) Wk{x)Z k j0 (x) must be a gradient of a spatial scalar. 
Therefore the function W k (z e ) in Eq. (7.17) must be the gradient of a function — @(z e ): 

W k {x-Z L ] = -@ >k {Z\x)) . (7.18) 

The condition (7.18) enables us to transform to zero the particular solution (7.17) of the 
Heisenberg equations of motion by successive applications of the symmetry transformation 
(3.27) generated by the Noether charge Q\&\. This transforms the momentum P k (x) = 
—P(x)W k (x) to zero. 

To summarize, any particular solution (7.18) of the constraint equations and the Hamil- 
ton equations of motion of non-rotating dust can be cast into the form in which the dust 
time foliation T± is orthogonal to the dust frame S. When we demand from the outset not 
only that the dust be non-rotating, but also that the time foliation T which is used for its 
description be orthogonal to the dust frame, we can impose the additional constraint 

P k (x) = 0. (7.19) 

This constraint is stronger than the requirement (7.10) that the vorticity fl a b(x) vanish; 
Eq. (7.19) implies Eq. (7.10), while Eq. (7.10) does not automatically enforce Eq. (7.19). 

Now consider the system of constraints H^(x) = 0, H^ k (x) = 0, and Pk{x) = 0. This con- 
straint system is first class: In Sec. 4 we showed that the new super-Hamiltonian constraints 
H^(x) and the new supermomentum constraints H^ k {x) have vanishing Poisson brackets 
among themselves. The momenta P k (x) also have vanishing Poisson brackets among them- 
selves and with the generators H^(x). From Eq. (4.9) we see that the Poisson brackets of 
P k and H^ k can be written as 

{P k (x), H n (x')} = Z1{x'){P k {x), H a (x')} + {P k (x), Zt{x')}H a {x') . (7.20) 

Because P k (x) are scalars under DiffS, their Poisson brackets with H a (x) are linear in P k (x). 
Thus, the Poisson brackets {P k (x), H^e(x')} also close. 

In general, the functional Q[&] from Eq. (3.25) does not have vanishing Poisson brackets 
with P k (x), and thus does not generate a symmetry of the theory defined by the new system 
of constraints = 0, H^ k = 0, and P k = 0. Indeed, for $ k ^ 0, Q[&] generates a tilt 
of the time foliation T that spoils the Gaussian character of the coordinate system (r, z k ) 
based on the irrotational dust. For -d = const, Q[&] generates a symmetry that consists of 
a relabeling of the leaves of the privileged foliation T±. In this case, Q[&] is a conserved 
charge whose value (for i9 = 1) is the total relativistic mass of the dust. From Eq. (3.24), 
the functional Q[$] is seen to equal the new constraint (7.19) smeared with a prescribed dust 
space vector field — £ k (Z(x)). Thus, the Poisson brackets of Q[£\ with P k (x) vanish weakly. 
The smeared constraint Q[$] is the canonical generator of dust space diffeomorphisms, in 
which the Lagrangian coordinate labels for the dust particles are changed. 

The additional constraint (7.19) enables us to replace the Hamiltonian constraint 
by an equivalent constraint that does not involve the square of the gravitational super- 
Hamiltonian H^_. Modulo Eq. (7.19), the new supermomentum constraint (4.9) implies 
H£ pa -\fGT ta . This can be combined with the definition G := (Hf ) 2 - g ab H^H^ to yield 
[Hff = (1 + g ab T a T b )G. The Hamiltonian constraint F T = P - \fe then takes the form 
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# T0 := P(x) + h (x; T, g,p] = , (7.21) 



with the Hamiltonian 

h (x; T,g,p] := (l + g ab (x)T a (x)T b (x)y 1/2 Hf(x) . (7.22) 

The new Hamiltonian (7.22) does not contain an inconvenient square root and is quadratic 
in the gravitational momentum p ab (x). Unlike the general Hamiltonian h(x) of Eq. (4.7), 
however, it does depend on the Gaussian time T(x) labeling the leaves of T±. The constraint 
(7.21) was employed for the quantization of gravity interacting with non-rotating dust in 
our previous work [10]. 



C. Quantizing a Non Rotating Dust 

We have seen that any particular solution of the classical constraint equations and the 
Hamilton equations of motion for non-rotating dust can be described with respect to the 
time foliation 71 that is orthogonal to the dust frame Z k ^ a (y). By imposing the constraint 
Pk{x) = 0, we demand both that the dust be non-rotating and that its motion be described 
with respect to the frame-orthogonal foliation. The canonical description of this system 
can be formulated either in terms of the original canonical variables or in terms of the new 
canonical dust-space variables. For the original variables, the new constraint has the simple 
form Pk{x) = 0, while the supermomentum constraint has a relatively complicated form 
(4.9). With the dust-space variables, the supermomentum constraint is simple — it is the 
momentum conjugate to Z h (x) — while the new constraint (7.19) takes the form 



P k (x) = H^ k (x) - 



dZ(x) 



dx 



H G k (z)+P(z)T^z] 



(7.23) 

--Z(x) 



In each case, the constraint that generates diffeomorphisms of the manifold on which the 
tensor fields live (S or S) has a relatively complicated form, while the constraint that 
generates diffeomorphisms of the complementary manifold is trivial. Thus, we see that for 
non-rotating dust with the frame-orthogonal foliation there is no relative advantage in using 
one set of canonical variables as opposed to the other. In what follows we will work with 
the original canonical variables. 

The system of constraints = 0, if-fo = 0, and P k = is first class. Therefore these 
constraints can be imposed consistently as restrictions on the state functional ^>[Z,T,g]. 
The constraint (7.19) in particular yields the restriction 

P k (x)*[Z,T,g] = 0, (7.24) 

and implies that the state \I/ does not depend on the frame variables Z k (x). Our task then is 
to find solutions ^f[T,g] of the operator equations that follow from the classical constraints 
H-\o(x) = and H^(x) = 0. These are the many-fingered time Schrodinger equation 

tf T o(x)*[T,(?]=0 ^ i-j^ = h i) (x ] T,g,p]*[T,g], (7.25) 
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and the equation that requires \& to be invariant under DiffX, 



T a (x)P(x) + H^(x))^[T,g] = . (7.26) 

Functionals ^f[T,g] that satisfy Eqs. (7.25) and (7.26) describe the quantum states of non- 
rotating dust with respect to the frame-orthogonal time foliation Tj_. We have explored this 



scheme in a previous paper 10 



One might expect that any state \l/ with vanishing vorticity, 

n ab (x)^[Z,T,g]=0 , (7.27) 

can be transformed into a state satisfying Eq. (7.24) by the action of the symmetry operator 
(6.25). However, this is not so. The underlying reason is that a classical Pk{x) can be 
brought to zero by a symmetry transformation for any particular solution of the equations 
describing a non-rotating dust, but the state functional *ff[Z, T, g] generally describes an 
ensemble of quantum systems corresponding to a superposition of such individual solutions. 
Correspondingly, there is in general no single function @(Z(x)) whose charge Q[&] would 
turn a state functional that satisfies Eq. (7.27) into one that satisfies the stronger equation 
(7.24). 

We now ask whether every Z-independent solution *S?[T, g] of Eqs. (6.47) is also a solution 
of Eqs. (7.25) and (7.26), and inversely, whether every solution of Eqs. (7.25) and (7.26) is 
also a solution of Eqs. (6.47). The proof of such a statement stumbles over a factor ordering 
problem. Assume that ^[T, g] solves Eq. (6.47a) and hence also the iterated equation 

P\x) *[T,g] = G(x; g,p] *[T,g] . (7.28) 

//the operator G(x) were factor ordered as 

G(x; g,p] = (H?(x)) 2 - g ab (x)H^(x)H^(x) , (7.29) 
then its application on a state ty[T,g] which also satisfies Eq. (6.47b) would yield 

(l +g ab (x)T a (x)T b (x))P 2 (x) * = (H^x)) 2 * . (7.30) 
The state \l/ would then satisfy Eqs. (7.25) and (7.26) with the operator h (x) factor ordered 



as 



h (x; T,g,p] = (l +g ab {x)T a {x)T h {x)) 1 Hf{x) . (7.31) 

However, the factor ordering of the operator G(x; g,p] is dictated by the requirement that 
the operators G(x) and G(x') commute, thereby replicating the vanishing of the correspond- 
ing classical Poisson brackets, Eq. (4.22). One can hardly expect that any ordering of the 
factors g and p in the individual operators Hl(x; g,p] and H^(x; g,p], followed by the or- 
dering (7.29) of G(x), can achieve this aim. Similarly, one cannot expect that the ordering 
(7.31) would yield Hamiltonian densities ho(x) that are self-adjoint. 

Formal considerations like these indicate that, even if we were able to construct con- 
sistent quantum field theories based respectively on the Schrodinger equations (6.47a) and 
(7.25), these theories would likely be non-equivalent. 
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APPENDIX A: NOTATION 

Manifolds. Manifolds play an important role in our study: spacetime At, space E, 
dust space S, and dust time T . Points y G At are labeled by local coordinates y a which 
carry Greek indices. Points x G E are labeled by local coordinates x a with indices from 
the beginning of the Latin alphabet. Points z G S are labeled by local coordinates z k with 
indices from the middle of the Latin alphabet. Instants r G T are always labeled by the 
Cartesian coordinate which (a minor ambiguity) is also called r. 

Mappings. Geometric properties of spacetime and physical properties of the dust are 
studied in the canonical formalism on (spacelike) embeddings Y :Y> ^ M.hy x EY* ^ y = 
Y(x) G Ai . The spacetime metric 7 on Ai is pulled back by the embedding mapping Y 
into the spatial metric g on E . The mapping Y : T x <S — > .M by (r G T, z <E S^J 1— > y = 
T(t, z) G Ai foliates Ai into the dust time T and the dust space S. The inverse mappings 
are T : M ^ T by y e M ^ t = T(y) GT and Z : M ^ S by y E M ^ z = Z(y) G S . 
A point x G E of the embedded hypersurface Y{Yi) lies at the instant r = T o Y(x) of the 
dust time and at the dust space point z = ZoY{x) . The proliferation of symbols is avoided 
by a slight misuse of notation: the mappings ToF:S^T and Z o Y : E — > S are simply 
called T and Z. We thus write T : E ^ Tby j; e S k r = T(x) G T and Z : E -> S by 
x G E 1— > z = G 5. It is fairly clear from the context what mappings the symbols T 

and Z denote. 

The mapping Z : E — > 5 has the inverse X = :5^Eby2;G5i-^x = -^(z) £ ^ • 
Tensor fields on E are mapped to S by Z*, and tensor fields on 5 are mapped to E by X*. 
The corresponding tensor fields are given lightface kernel symbols on E and boldface kernel 
symbols on S. For example, the spatial metric g a b(x) on E is mapped to g k i(z) on S and 
the dust time T(x) on E is mapped to T(z) on S . 

Notice that in all these examples, the mappings are denoted by capital letters, and their 
values by the corresponding lower case letters. 

Derivatives. Partial derivatives are denoted by comma followed by the appropriate 
coordinate index: >a , jfl and ^ . The covariant derivative on (Ai,^) is written as V a , the 
covariant derivative on (S,g) as D a . The pound symbol £ is used for the Lie derivative. 

Brackets. The arguments of functions are enclosed in round brackets ( ), the arguments 
of functionals in square brackets [ ] . Mixed brackets ( ] indicate that an object is a function 
of one argument and a functional of another: e.g., R(x; g] means that the scalar curvature 
R on E is a function of x G E and a functional of the induced metric g a b(x') ■ 

The Lie brackets of vector fields and the commutators of operators are denoted by square 
brackets [ , ] . Braces { , } are reserved for the Poisson brackets. 

Densities and operators. When the same quantity appears both as a tensor and as 

a tensor density, the density is marked by a tilde. Thus, J a is a vector while J a is a vector 
density of weight 1. All operators wear hats, as g(x) or p(x). 
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Units. We are using natural units in which the speed of light c and the rationalized 
Newton's constant 16nG are set equal to 1 : c = 1 = IQ-kG . In these units, the gravitational 
Lagrangian is simply the densitized scalar curvature |7| 1//2 i?(y; 7], and the Einstein law of 
gravitation takes the form R a/3 - \Rf0 = \T a ^ . 

APPENDIX B: ALTERNATIVE FORMS OF THE SPACETIME ACTION 

The action (2.4) is a functional of eight state variables, T, Z k and M, W k . It can be cast 
into an equivalent form that depends on fewer state variables by solving the Euler equation 
(2.7), obtained by varying W k , with respect to W k . This yields 

W h = huZ'^Tp , (B.l) 

where h k e is the inverse of the matrix 

h u := Z\ al a ^Z\p . (B.2) 

The metric hu measures the distance ds orthogonal to the flow lines between neighboring 
flow lines with Lagrangian coordinates z k and z k + dz k : 

ds 2 = h M dz k dz e . (B.3) 

By using Eq. (B.l) to eliminate W k from the action, we obtain the equivalent action 

S D [T, Z k ; M; laP \ := S D [t, Z k ; M, W k = W k [T, Z\ 7^]; 7q/3 ] (B.4) 

which depends on the four dynamical variables T, Z k , and a single multiplier M. 

The mass multiplier cannot be eliminated at the spacetime level by an algebraic process, 
because the equation (2.6) obtained by varying M does not depend on M. In Sec. 3a 
we showed how M can be eliminated from the canonical action (or, rather, how M can 
be replaced by the canonical momentum conjugate to the time function T). Instead of 
decreasing the number of variables by solving the field equations, one can increase the number 
of variables by adjoining valid equations to the action by additional Lagrange multipliers. 
As an example, the definition (2.2) of U a can be adjoined to the action (2.4) by a set of four 

multipliers J a that are densities of weight I. One obtains thereby the action 

S D [T, Z k ;M,W k ;U a ,J a ; la p] 

= J m d 4 y ^-i|7| 1/2 M(7 a/3 f/ a ^ + 1) + J a (U a + T a - W k Z k ^ (B.5) 

which depends on sixteen state variables T, Z , M, W k and U a , J a . By varying the action 
(B.5) with respect to J a , we recover the definition (2.2) of U a . By varying U a , we obtain 

J a - \-f\ 1/2 MU a = . (B.6) 
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This shows that the multipliers J a are the components of the (densitized) mass current 
(2.15). By calculating U a from Eq. (B.6), 

U a = M- 1 \^ 2 J a , (B.7) 
and substituting it back into the action (B.5), we obtain an equivalent action 

S D [T, Z k ;M,W k ;J a ; la p] 

= J M d 4 y ^r^M'^rj? - \\i\ l,2 M+ J a (T a - W k Z\ a )\ (B.8) 

that depends on twelve state variables T, Z , M, W k and J a . The variation of the action 
(B.8) with respect to M yields an equation that determines M, namely, 



M=\ 1 \- 1/2 \j- la pJ a jP . (B.9) 
By using Eq. (B.9), M can be eliminated from the action (B.8). The result is the action 

S D [T, Z k -W k J a - la(3 \ 



M 



d A y -y-iapJ a J (3 + J a (T a - W k Z\ a ) (B.IO) 



that depends on eleven state variables T, Z , W k and J a . This action for dust was introduced 



in Ref [26] in the study of perfect fluid actions. 

All of the actions (2.4), (B.4), (B.5), (B.8), and (B.IO) lead to the same canonical action. 
We have chosen the action (2.4) because the number of state variables exactly matches the 
number of canonical coordinates plus conjugate momenta. The physical meaning of the 
momenta then clearly emerges from the Legendre transformation. 



APPENDIX C: SYMMETRY TRANSFORMATIONS 

Transformations of the dust variables (T, Z k , W k ) — > (T", Z k \ W k >) that satisfy 

S D [T'(T, Z l , W e ), Z k '(T, Z\ W t ), W k ,(T, Z e , W e )] = S D [T, Z k , W k ] (C.l) 

constitute symmetries of the dust action. Inspection of the action (2.4)-(2.5) shows that 
transformations that preserve the Pfaff form (2.3), 

- dT + W k dZ k = -dT' + W k ,dZ k ' , (C.2) 

are symmetries. A Pfaff form —dT + W k dZ k , with k = 1, 2, 3, always can be reduced to 



dA + B dC for some functions A, B, and C ||18[|; however, here we are interested in invertible 
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transformations that leave the action functional invariant, as expressed though Eq. (C.l). 



(See Refs. [26,27 for discussions of the use of the reduced Pfaff form in the construction of 
perfect fluid action functionals.) We shall confine our attention to transformations that are 
ultralocal in spacetime, so the primed variables at a given spacetime point y depend only on 
the unprimed variables at the point y. In particular, the primed variables are not allowed 
to depend on derivatives of the unprimed variables. 

For the case of ultralocal transformations, the relationship (C.2) can be written as 

W V {T, Z\ W e ) dZ k '(T, Z e , W t ) = W k dZ k + dT{T, Z l , W e ) (C.3) 

and viewed as a time T-dependent canonical transformation for a fictitious phase space 
with coordinates Z k and momenta W k PS|,Z7||. Here, T is the function of T and of the 



'canonical variables' Z k , W k , defined by JF := T' — T. Now, for those transformations 
that are connected to the identity, the dependence Z k = Z k (T, Z , We) of the primed 
'coordinates' on the unprimed variables can be inverted for the unprimed 'coordinates', 
Z k = z k (T,Z e ',W e ). Then T, Z k ' , and W k can be treated as independent variables by 
substituting the functional relationship Z k = Z k (T,Z e ,Wi) into Eq. (C.3). This leads to 
the expression 

W k ,(T, Z 1 ', W t ) dZ k ' = W k dZ k (T, Z e ',W e ) + dF(T, Z l \ W e ) , (C.4) 

with a slight abuse of notation in which W k /(T, Z e (T, Z m \W m ), We) is set equal to 
W k '(T, Z e , We), and likewise for T . In terms of the generating function F(T, Z e , We) := 
JF(T, Z e , W e ) + W k Z k (T, Z* , We), Eq. (C.4) can be rewritten as 

W k '(T, Z e ', We) dZ k ' + Z k {T, Z e ' ,W e ) dW k = dF(T, Z e ' , We) . (C.5) 

This expression along with the definitions of T and F yield the results 

T '- T=F - W *iw i < (a6a > 
z " = lk< (a6b) 



where F = F(Z e ' , We) is independent of T. 

The Noether currents (2.29) associated with one-parameter subgroups of symmetries 
F{Z 1 ' , We, a) connected to the identity F{Z 1 ' , We, 0) = Z k 'W k are given by 

J a [f] = \l\ 1/2 J a f(Z e ,W e ) , (C.7) 
where J a is the mass current (2.15) and 

dF(Z e ',We,a 



f(Z\We). 

oa 



(C.8) 

CT=0 
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The corresponding Noether charges (2.28) can be written as 

Q[f] = -J :s d 3 x\g\ 1 ' 2 n a J a f(Z e ,W e ) (C.9a) 
= J^d 3 xPf(Z e ,-P £ /P) . (C.9b) 

In expression (C.9b), we have made use of the canonical variables P and P k as defined by 
Eqs. (3.10) and (3.12). 

In Sec. 2e of the main text we restrict our attention to those symmetry transformations 
that mix the dust spacetime variables T and Z k with one another. That is, we consider 
transformations for which T' and Z k are functions of T and Z k only. If this restriction is 
imposed on the transformations (C.6), it follows that the most general generating function 
F has the form 

F(Z e , W e ) = W k ~ k (Z e ') + @(S e (Z m ')) . (CIO) 
The corresponding symmetry transformations (C.6) are given by 

T' = T + @{Z m ), (C.lla) 

Z k = S fc( Z m') ^ (C.llb) 

W v = E\ k ,(Z m ')(w e + &AZ m )) , (Cllc) 

which coincide with the transformations (2.26)-(2.27) discussed in Sec. 2e. 

The Lie algebra of one-parameter groups of symmetries generated by F of the form 
(C.10) is labeled by the functions 



H . dF(Z e ',W e ,a) 



f(Z',W e ) := 



da 



= W k e(Z e ) + #(Z e ) , (C.12) 



cr=0 



where £ k (Z e ) and iD(Z e ) are defined in Eq. (2.32). Insertion of this expression for / into the 
Noether charge (C.9) yields 

Q[£,0] = ~jj^ \g\ 1/2 n a (j%e(Z e ) + J a #(Z e )) (C.13a) 
= Jj^x (-P k e(Z £ ) + P^(Z e )) , (C.13b) 

where J a and J k are the mass current (2.15) and the momentum current (2.16), respectively. 
The charges <3[£,#] coincide with the Noether charges (2.36)-(2.37) and (3.24)-(3.25) pre- 
sented in the main text. 



APPENDIX D: THE ACTION FOR RELATIVISTIC PARTICLES 

The geodesic motion of a single relativistic particle of rest mass m is generated by the 
super-Hamiltonian 
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H:=±-(^{y)P a P p + m 2 ) . (D.l) 

The spacetime coordinates y a of the particle and the four-momentum P a are canonically 
conjugate variables. These variables can be split into their space and time components with 
respect to the foliation (3.1)-(3.2). The momenta conjugate to t and x a are 

Pt :=X a P a and p a :=X a , a P a . (D.2) 

By using Eqs. (3.3) and (3.4), the super-Hamiltonian constraint H = can be resolved with 
respect to the energy variable — p t : 

~Pt = tl(t,X a ,p a ) 



-N a {t, x)p a + N^t, x)Jm 2 + g ab {t, x)p a p b . (D.3) 



The expression h(t,x a ,p a ) is the Hamiltonian of the particle. 
The action associated with the Hamiltonian (D.3) is 

S[x a , Pa ] = jdt L(x a ,p a ,x a ) = Jdt (p a ^- - hit.x^pa^j . (D.4) 

Observe that the derivative of the canonical Lagrangian L with respect to mass m is 

N± . (D.5) 

dm yj 1 + g ab p a Pb/m 2 

Since N ± dt is the orthogonal proper time separation between t = const surfaces, and 



1 + g ab p a Pb/™ 2 is the relativistic 'gamma factor', dL/dm equals the ratio dr/dt, where r 
is the proper time along the particle worldline. This shows that the action (D.4) can be 
extended to a functional of m(t) and r(t), 

S[r,m; x a ,p a ] = J dt(m^ + p a ^ - N a (t, x){-p a ) 

-N^t, x)yjm 2 + g ab (t, x)p aPb ) . (D.6) 

Extremization of S with respect to m relates dL/dm to dr/dt. Extremization of S with 
respect to r yields the condition dm/dt = 0, that the particle mass is a constant of motion. 

We now compare the relativistic particle action (D.6) with the canonical dust action 
(3.16). For the purpose of comparison, the metric tensor components N ± (t,x), N a (t,x), 
g a b{t,x) are treated as prescribed external fields. From the analysis of Sec. 5 it follows that 
the dust action (3.16) can be written as a functional of canonical variables T(z), P(z), 
Z k (x), and H^ k (x) := Z%(x)H® (x). We can make a further canonical transformation [15] by 
replacing the dust frame Z k (x) with its inverse, X a (z). The new variable X a (z) specifies the 
spatial location of the dust particle with Lagrangian coordinate label z k . The momentum 
conjugate to X a (z) is found from the relationship 
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Jj 3 xH^ k (x)Z k (x) = -Jj 3 z H^ k {X{z))Z\ a (x{z))x a {z) . (D.7) 

With the momentum conjugate to X a (z) denoted by P a {z), the action (3.16) becomes 
S D [X a ,P a ; T,P]= f dt f d 3 z(p(z)f(z) + P a (z)X a (z) - N a (t,X(z)) (-P a (z)) 

J IR J S \ 



-N L (t,X(z))^J 'p\z) + g<*(t,X(z))P a (z)P b (z) . (D.8) 



This action coincides with the relativistic particle action in the following sense. The action 
for many relativistic particles is obtained from the single particle action (D.6) by summing 
over the individual particles. With the particles labeled by z = 1,2,..., the canonical 
variables are given by the set (r 2 , m z , x a z) p za )- In the continuum limit of an infinite number 
of particles scattered throughout space, the sum J2 Z over individual particles is replaced 
by an integral / d 3 z over the space of particle labels, the dust space S. At the same time, 
the canonical momenta m z and p za become densities in the space of labels. The resulting 
action functional is precisely the dust action (D.8), with the correspondences r z — > T(z), 
m z — > P(z) d 3 z, x a z — > X a (z), and p za — > P a {z) d 3 z. 



APPENDIX E: THE POISSON ALGEBRA OF G(x) 

We check by direct calculation that the quadratic combinations G(x), Eq. (4.7b), of 
the gravitational super-Hamiltonian if^(x) and supermomentum H^(x) have strongly van- 
ishing Poisson brackets, Eq. (4.22). The result follows from the fact that the gravitational 
generators H^_(x) and H^(x) satisfy the Dirac 'algebra' Eqs. (4.1)-(4.3), and from the action 
of these generators on the metric g ab (x). 

Let us first show that the expressions 

F(x) := g ab (x)H^x)H^(x) (E.l) 

commute among themselves, 

{F(x),F(x')} = . (E.2) 
From the definition (E.l), it holds that 

{F(x),F(x')} = (2{F(x),H°(x')}H Ga (x') - (x <-> x')) 

-4{HC(x), H^x')}H g \x)H g \x') . (E.3) 

Now smear Eq. (E.3) by two externally prescribed scalar densities of weight —1, M(x) and 
N(x'). With the abbreviations 

M a (x) := M(x)H Ga (x) , N a (x) := N(x)H Ga (x) , (E.4) 

Eq. (E.3) becomes 



47 



{F[M],F[N]} = (2{F[M], H G [N]} - (N <-> M)) - 4{# G [M], H G [N]} . (E.5) 

Note that in Eq. (E.5), M and N must be treated as externally prescribed smearing vectors. 
This is because the supermomenta H Ga (x), which combined with M(x) and N(x) into the 
smearing vectors (E.4), stand outside the Poisson brackets in Eq. (E.3). 

Since F(x) is a scalar density of weight 2 under DiffE, its Poisson bracket with H G [N] 
is given by 

{F[M],H G [N}} = J^d 3 xM(x)£^F(x) 

= f d 3 xM (F a iV a + 2FN a , a ) . (E.6) 

Therefore the first term on the right-hand side of Eq. (E.5) yields 

2{F[M],H G [N]} - (N <-> M) = aJ d 3 x F (MN a - NM, a ) H Ga . (E.7) 

The second term on the right-hand side of Eq. (E.5) can be evaluated from the smeared 
form of Eq. (4.3), which asserts that H G [M] represent the Lie algebra DiffS: 

{H G [M],H G [N]} = H G [[M, N}] . (E.8) 

Here, [M, N] is the commutator of the vector fields M and N, 

[M, N} b = N\ a M a - M\ a N a 

= (MN A - NM A )H Ga H Gb . (E.9) 

Equations (E.8) and (E.9) imply that the second term on the right-hand side of Eq. (E.5) 
exactly compensates the first term given by the integral (E.7). This proves Eq. (E.2). 
The Poisson bracket that interests us is 

{G(x),G(x>)} = {F(x),F(x>)} 

-(2H G (x){H G (x),F(x')}-(x~x')) 
+AH G (x)H G (x'){H G (x),H G (x')} . (E.10) 

We have just proved that the first term on the right-hand side of Eq. (E.10) vanishes. 
Because H G (x) induces an ultralocal change of the metric, 

{g ab (x), H G (x')} = 2K ab (x)5(x, x') (E.ll) 

(the coefficient K ab is the extrinsic curvature), the action of H G on g ab in F cancels out 
under the interchange x <-> x'. Therefore, by Eq. (4.2), the second term on the right-hand 
side of Eq. (E.10) yields 

-AH G (x){H G (x),H G (x')}H Ga (x') - (x <-> x') 

= AHf(x)Hf(x')H Ga (x')5 tal {x, x') - (x <-> x') . (E.12) 

Up to a sign, this is what the last line of Eq. (E.10) yields when we substitute into it the 
Poisson bracket (4.1). As a result, the right-hand side of Eq. (E.10) vanishes. We have 
thereby checked the pivotal equation (4.22) of the main text. 
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APPENDIX F: DiffiS AND THE EIGENFUNCTIONALS OF G{z) 
We want to prove that if <p[g] is an eigenfunctional of G(z) with eigenvalue G(z), 

G(z) <t>[g] = G{z) 4>[g] , (F.l) 
then 0[3*<jr] is an eigenfunctional of G(z) with eigenvalue (H _1 )*G(,z) : 

G(z) <j>[S*g] = (S-'TGiz) cj>l~*g] . (F.2) 
It is sufficient to prove the infinitesimal version of Eq. (F.2), namely 



d_ 

da 



G{z)<t>[K9. 



cr=0 



d_ 

da 



"z- l yG{z)<t>[Kg\ 



CT=0 



(F.3) 



and then to pass to an arbitrary finite diffeomorphism H by exponentiation in a. 

To begin, recall that the classical charge Q[$] discussed in Sec. 5b is the canonical gen- 
erator of DifLS. Thus, the action of DiffiS on the operator G(z) is given by the commutator 



-[G(z),Q[^]] = £ i G(z) . 



(FA) 



Also observe that the change in <p[g] under an infinitesimal diffeomorphism of its argument 
g is determined by the operator Q[£] acting on (f>[g]: 



<T=0 



= i I d?z(£- i g u (z))p u (z)^[g\ 

J s 



(F.5) 



By letting G(z) act on equation (F.5) and using the commutator (F.4), we obtain 
d 



Ta [G ( z)m'A 



CT=0 



= i[G(z),Q[£]](f>[g]+iQ[Z]G(z)<f>[g] 



(F.6) 



Now use Eq. (F.5) along with the fact that <j>[g] is an eigenfunctional of G(z). Equation 
(F.6) then becomes 



d 

da 



G(z)<p[s* a g} 



CT=0 



-f^-rG(z) 



da 



^_4>[g] + G(z)^-[4>[K9) 



<T=0 



;f.7) 



which is equivalent to Eq. (F.3). This proves the theorem expressed in Eq. (F.2) 
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APPENDIX G: RELATIVISTIC PARTICLE IN THE PROPER TIME 

FORMALISM 

As in Appendix D, let us regard the rest mass m for the relativistic particle as a dynamical 
variable canonically conjugate to the proper time r along the particle's worldline: 

{r,m} = l. (G.l) 

The proper time formalism for the relativistic particle is obtained by resolving the super- 
Hamiltonian constraint in the extended phase space (r,y a ; m, P a ^j , 

H cxt := m 2 + 1 aP {y)P a P p = , (G.2) 
with respect to m. This yields 

/ lT ;= m + / l = 0, h := - s l- 1 cP(y)P a P p , (G.3) 

where the choice of the negative square root in the Hamiltonian h ensures the positivity of 
m. The reality of m requires that 

G := -r P (y)P a Pp > , (G.4) 

i.e., that the particle moves along a timelike worldline. 

The imposition of the constraint (G.3) as an operator restriction on the state function 
\I/(t, y a ) yields the Schrodinger equation in proper time r: 

H^ = ^ id T *(T,y) = h*(T,y) . (G.5) 

The Hamiltonian 

h = -y/^6, 6 = ^ (G.6) 

must be defined by spectral analysis. It is self-adjoint on the Hilbert space Tt + spanned by 
the eigenfunctions (fie of G belonging to positive eigenvalues G : 

GMy) = GMy) , g>o. (g.7) 

On TC + , the Schrodinger inner product 

:= J d 4 y^(T,y)^(r,y) (G.8) 

is conserved in r. Observables must leave the states in H + , i.e., they must satisfy the 
condition (6.46), where 

V + = 0(G) (G.9) 

is the projector to H + . 
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In general, neither y a = y a x nor P a = —id/dy a are observables. This prevents us 
from interpreting the multiplication operator y a as an observable whose measurement at the 
proper time r would localize the particle at a spacetime event y a , and ^(r, y)^(r, y) as the 
probability density that the localization will occur at y a . 

This difficulty persists even for a particle moving in flat spacetime ^^{y) = f] a ^ ■ In 
this case, P a commutes with G and hence it is an observable. On the other hand, one can 
explicitly show that y a = y a x is not an observable since one can always choose a constant 
covector 7c a such that the action of the unitary operator 

U w :=exp(in a y a ) (G.10) 

leads the state \I> out of 7i + . 

The demonstration follows the general outline of an old argument by Pauli [53 that 
if time were represented by the multiplication operator, the energy spectrum could not 
be bounded from below. In our example, the Hilbert space 7i is spanned by the joint 
eigenfunctions 

<S> P (y) = (27i)- 2 exp(tP a y a ) (G.ll) 
of the momentum operators P a . Equation 

*(j/) = Jd 4 P^(P)$ P (y) (G.12) 

relates the ^-representation $?(y a ) of a state \I/ G 7i to its P-representation ^f(P a ) . The 
momentum eigenstates (G.ll) are also eigenstates of G with eigenvalues G = —rj a ^P a Pp. 
The state \l/ lies in 7i + if the state function ^(P a ) has its support only on timelike vectors 

*(P a ) = for f] a ^P a Pp > . (G.13) 

The action of U n on the state function ^(P a ) translates its argument by —7i a : 

(U^)(P a )=^(P a -ir a ). (G.14) 

The support of (U^) (P a ) is thereby shifted by the amount n a relative to the support of 
^(P a ). Because ir a is arbitrary, the support of (U n ^) (P a ) cannot always remain timelike. 
This means that y a is not an observable. 

This simple example illustrates why the fundamental gravitational variables cannot be 
expected to be observables. 
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Z k (y) = z k = const Z k (y) = z k + dz k = const 
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T(y) = x + dx = const 
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v v " 
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Figure 1: U a d/dy a = d/dr and T a tk d/dy a = d/dz k are the coordinate basis vectors for 
the dust coordinates (r, z k ). The covectors T^ a dy a = dr and Z k ^ a dy a = dz k correspond 
to the stacks of surfaces r = const and z k = const, respectively; the associated vectors are 
T ,a ._ g af3 T/3 and Z k,a ._ g ^z k ^. The vectors U a and Z k > a are orthogonal, U a Z k , a = 0, 
and are therefore drawn symmetrically about the light cone. Likewise, the vectors T' a and 
T a ik are orthogonal, T Q T a ifc = 0, and are drawn symmetrically about the light cone. 
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Figure 2: The hypersurface Y (t, S) is displaced to the hypersurface Y(t + dt, E) along the 
dust flow lines z k = const by the proper time dr — N^(x — Z^ 1 (z)^jdt . The displacement 
is generated by the Poisson brackets with H[N^] dt. 
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Figure 3: The hypersurface Y(t, S) is displaced to the hypersurface Y(t + dt, X) along the 
leaves T(r, S) of the dust time foliation T by the amounts dz k = N^^x = Z~ 1 (z)^J dt in the 

dust space S. The displacement is generated by the Poisson brackets with H[N^} dt . 
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